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ELECTRON AND ION TEMPERATURES IN THE IONOSPHERE* 

A. P. Willmoret 

Physical Research Laboratory, Ahmedabad, India 

and 

Laboratory for Planetary Atmospheres, Goddard Space Flight Center 

Greenbelt, Maryland 

ABSTRACT 

Many measurements of ionospheric temperature have been made during the 
last decade and have shown a considerable departure from thermal equilibrium 
there. Theoretical work has led to a general understanding of the processes 
which determine the energy balance of the charged particles , and most features 
of the experimental results are well understood in the F -region. The position 
in the E-region where different methods of measurement, known to agree well 
at greater heights, appear to disagree is less clear, whilst in the D-region no 
direct observations of charged particle temperature exist. The theoretical ex- 
pectation here is that, on account of the rapid exchange of energy amongst the 
particles by collisional processes, temperature equality will prevail. 

*This review is based on material published up to January, 1970. 

t On leave of absence from Mullard Space Science Laboratory, University College, London. 
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ELECTRON AND ION TEMPERATURES IN THE IONOSPHERE 


1. INTRODUCTION 

The first direct attempt at measurement of electron temperature in the 
ionosphere was made in 1947 by Reifman and Dow (1949). Subsequently, a large 
number of results has become available from rocket, satellite and ground-based 
observations, and it is found that generally the charged particle temperatures 
exceed, often by a considerable factor, that of the neutral particles. There are 
three main reasons why these results are of interest. The charged particle 
temperatures are important parameters for the distribution of ionisation in the 
topside ionosphere, where diffusive equilibrium prevails, and influence the rate 
of recombination at lower altitudes. The redistribution of ionisation following 
temperatime changes may contribute in an apparently anomalous way to the di- 
urnal changes at the layer peak, and changes in recombination rate during iono- 
spheric storms may also be important. The major source of energy which main- 
tains the temperature difference with the neutral gas, with which the charged 
particles are in thermal contact as a result of collisions, is solar EUV radiation, 
but the possibility exists of revealing the existence of other sources at night, at 
high latitudes or during ionospheric storms. Finally, the close relationship be- 
tween the magnetosphere and ionosphere has to a considerable degree been 
recognised as a result of the study of charged particle temperatures and the 
diurnal interchange of energy between these regions . 

After the pioneering but unsuccessful attempt of Reifman and Dow to apply 
in the ionosphere the laboratory technique of the Langmuir probe, development 
continued in several laboratories resulting in successful results being obtained 
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in increasing quantity, both from rockets and from satellites, from 1959 onwards. 
During the same period the increasing availability of electron number density 
protiles for the topside ionosphere, obtained from topside sounding or by other 
techniques, lead to many analyses of these profiles using the diffusive equilib- 
rium theory and the deduction of plasma temperature from them. Soon after, the 
suggestion originally made by W. E. Gordon in 1958 that electron temperature 
might be measured by observing the Thomson scatter of a radio wave from the 
ionosphere lead to a rapid development of the incoherent scatter technique and 
the establishment at several locations in Europe and America of observing 
stations. 

This rapid development of experimental methods and the early discovery 
that at altitudes near 250 km at least the electron temperature appeared to be 
substantially above the neutral gas temperature stimulated an interest in the 
collisional processes which determine the charged particle temperatures. The 
original basis of the theory was laid by Hanson and Johnson in 1961 who showed 
that solar radiation, in ionising the neutral gas, deposited a fraction of its energy 
in the ionospheric electrons thus raising their temperature significantly. This 
theory underwent substantial development in the hands of several other authors 
and it is now clear that it accounts for the major features of the quiet F-region 
during the day and, surprisingly, is even able to account for a departure from 
thermal equilibrium which is observed to occur at night when solar radiation is 
absent. Interest now centres largely around the complicated changes which oc- 
cur in ionospheric storms, and effects which may be observed over the auroral 
zone. 
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In the D and E-regions the situation is unclear. Several direct comparisons 
of temperatures measured by different methods have been made, in particular 
between Langmuir probes and incoherent scatter, with excellent results above 
130 km. Below this altitude, no such comparison has been made and here rather 
a strange inconsistency between the results obtained by the two methods is found. 
It is thus not clear whether a departure from thermal equilibrium occurs in the 
E-region; if so, it cannot be sustained by the same collision processes as in the 
F-region as sufficient energy is not available. On the other hand, it is by no 
means impossible that solar radiation is the original energy source. At present 
no direct observations are available in the D-region. 

In this paper the principal observing techniques will first be briefly reviewed. 
Then the theory of solar heating of the ionosphere will be described and' com- 
pared with observations of the quiet ionosphere in the F-region. Then the situ- 
ation in the E-region will be described. Next observations of temperature vari- 
ability over the solar cycle and during ionospheric storms will be discussed, 
followed by an account of processes of special importance at high latitudes. 

2. TEMPERATURE MEASUREMENTS BY LANGMUIR PROBES 

Electrostatic probes were first used in measurements of the potential 
distribution in gas discharges by J. J. Thomson, their theory of operation later 
receiving extensive development through the work of I. Langmuir and his col- 
laborators [Langmuir (1923) (a), (b), (c), Langmuir and Mott-Smith (1924), 
Mott-Smith and Langmuir (1926)] . This led to extensive application of probes 
in the study of gas discharges, and, together with some later developments, 
paved the way for the application in the ionosphere where conditions are actually 
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more favourable for the utilisation of the method than m a gas discharge and 
where the most elegant applications of it have occurred. Although we shall be 
concerned with modes of operation of the electrostatic probe which Langmuir 
had not envisaged, it nonetheless seems appropriate to recognise his leading 
contributions to the theory by the use of the general term, Langmuir probe. 

The Langmuir probe is then an electrode of arbitrary geometry immersed 
in a charged particle plasma, whose current -voltage characteristic is determined. 
Prom it, under suitable conditions, the number densities of the ions and electrons 
and their energy distribution functions may be obtained. The current to the probe 
will an general arise from the collection of positive ions , negative ions and elec- 
trons from the plasma. Negative ions, however, are present in significant con- 
centrations only in the D-region where collisions near the probe due to the rela- 
tively high gas densities render the theory as discussed here inapplicable, so 
that the presence of negative ions will be ignored. Very few efforts have been 
made to produce a probe theory applicable in the D-region [Hoult (1965), Hoult 
and Ta-Jin Kuo (1966), Sonin (1967)] ; the results suggest that information on the 
energy distribution function or temperature cannot be obtained when collisions 
are important. Besides the collection of charged particles from the plasma, 
contributions to the probe current may also arise from the emission of electrons 
from the probe surface due principally to irradiation by sunlight and from the 
collection of secondary or photo -electrons from the spacecraft surface. 

The magnitude of each of these components of the current depends on the 
potential V p of the probe relative to that of the undisturbed ambient plasma. 

The basic problem of probe theory is to relate the observed current-voltage 
characteristic to the energy distribution function of the particles. The properties 
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of the regime in which the probe accelerates a given particle type are funda- 
mentally different from that where the particles are retarded, especially as far 
as the effect of the details of the electric field distribution about the probe are 
concerned. It is found that relations of more general validity exist for the re- 
tarding regime. In the accelerating regime, the results may be affected by the 
formation of a sheath about the probe. Immediately adjacent to the probe, the 
plasma will be deficient in particles for which V p is repulsive, leading to a 
region of net space charge density. This will also be a region of large electric 
field, so that in the presence of a plasma, the electric field strength is stronger 
close to the probe surface and weaker at large distances, than if the plasma is 
absent. The region of strong field is called the sheath; its effect is to screen the - 

plasma from the applied potential on the probe. The sheath thickness depends 

1 _ 

on V p but is usually (5 - 20)A. D , where A. D = (kT e /4-7rn e e 2 ) 2 is the Debye length 
in the ambient plasma. Obviously, the retarding and accelerating regimes occur 
simultaneously for the two types of plasma particle, and some experimental 
method of distinguishing their contributions as well as those of secondary or 
photoelectrons must be found. 

In practice the electron current is usually very much larger than the other 
contributions except high in the magnetosphere so that a fairly simple correction 
is sufficient. This can be made by determining the probe current at a large nega- 
tive potential where the electron current is effectively absent and using this value 
as the current zero. A more accurate determination, essential in the case of the 
ions, can be made by using a screened probe [Boyd (1968), Bordeau and Donley 
(1964)] in which grids are used to separate the different contributions to the 
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current. This system requires careful design however because the effect of the 
fringe fields near the grid wires can be much larger than might have been 
anticipated. 

In practice the Langmuir probe characteristic can be calculated only in 
idealised cases and there is m consequence a large number of potential sources 
of error. Intercomparison of probe measurements and of probe and incoherent 
scatter results, described in Section 6, suggests that these are mostly under- 
stood. In this section, we first derive the probe characteristic for certain cases 
and then discuss briefly the sources of error. 

(a) The Retarding Regime 

Consider for definiteness the electron current to a probe of arbitrary 
geometry immersed in a plasma. We shall suppose the electrons to have a 
Maxwellian velocity distribution in a coordinate system fixed with respect to 
the probe. It will be shown in Section 7 that in the ionosphere the electron ve- 
locity distribution is quite closely Maxwellian referred to a coordinate system 
which is co-rotating with respect to the Earth. The mean thermal velocity of 
the electrons exceeds by an order of magnitude even satellite velocities and 
Ragab (1968) has shown that the effect of spacecraft movement on the electron 
probe characteristic is not large. In the retarding rdgime, the electron density 
near the probe surface will be given by Boltzman's relation 

n P = n o exp(-eV p /kT e ) (1) 

where n p and n o are the electron densities near the probe surface and in the 
ambient plasma, respectively, and 
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T e is the electron temperature, 
e is the electronic charge and 
k is Boltzman's constant. 


This, of course, is an equilibrium relationship which will be applicable 
accurately only if the probe surface is perfectly reflecting. However, it will' be 
a good approximation provided the probe surface current density j is suffi- 
ciently small compared with the random current j o in the plasma. 


<< 
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where c e is the mean thermal velocity of the electrons. In fact, even with this 
restriction, there are important limitations on the applicability of Eq. (1) which 
are further discussed in Section 3(b) and (c). From (1) and (2) we then obtain 
Langmuir's expression for the electron retardation characteristic 

i p = Jo e*p(-eV p /kT.) (3) 


which may- be used to derive both n o and T e from a measured characteristic. 

For the ions, too, the assumption of a Maxwellian velocity distribution in 
co-rotating coordinates may be justified [though see Section 7 (a)] , but in the 
reference frame of the measurements the probe velocity must now be taken into 
account as it is always comparable with the mean thermal speed of the ions and 
may exceed it by nearly an order of magnitude. In this case, we may use a 
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generalisation of the above approach, due to Druyvesteyn (1930). If the velocity 
distribution function of the particles is isotropic, the ordinate of the energy 
distribution function being given by f (V), then the equivalent of Eq. (3) is 



3 2 V 2 

P 



( 4 ) 


which reduces to (3) in the case of a Maxwellian distribution. In our case, where 
the velocity distribution is anisotropic it is easy to show that this relation still 
holds for a spherical probe. 

Thus, from the probe characteristic in the retarding regime it is possible to 
find the energy distribution function (and hence, the number density) of the elec- 
trons or the ions. The results are quite independent either of the probe geometry 
or the precise distribution of the electric field about the probe and so of sheath 
formation, provided that the particle velocities are isotropic. Otherwise they 
are valid only for a spherical probe but are still unaffected by sheath formation 
if the electric field distribution remains spherically symmetric. This condition 
cannot be precisely satisfied for an anisotropic velocity distribution but appears 
to obtain with sufficient accuracy in practice. 

(b) The Attracting Regime 

When the probe potential is such that the particles are attracted to it, de- 
termination of the probe current requires that the particle trajectories them- 
selves be determined. This requires a detailed knowledge of the electric field 
distribution and hence of the sheath formation. No simple solution of the general 
problem is now possible. It will be instructive to consider the case of a spherical 
or cylindrical probe of radius a. if the Debye length \ >> a, the effect of space 
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charge can be ignored and the probe current i p is then found to be [Mott-Smith 
and Langmuir (1926)] 



where i D is the probe current when V p - o 
1 

and a =~ for a cylinder and 1 for a sphere 

V depends on the particle energy, being kT for a Maxwellian distribution and V o 
for a homogeneous stream of energy V o . Under these circumstances the probe 
current is orbital motion limited. 

If on the other hand A. D << a, the probe is surrounded by a thin sheath. The 
current density crossing the sheath edge will be j while because of the strong 
radial electric field in the sheath the particles entering the sheath will have a 
high probability of being swept to the probe. Hence, for all probe potentials 

l = i (6) 

and the probe current is sheath limited. 

For intermediate values of \ D , which is the practical case, the problem 
arises of what relation should be used in place of (5) or (6). This question was 
discussed by Mott-Smith and Langmuir (1926) and by Allen, Boyd and Reynolds 
(1957). In fact (5) is derived from rather general considerations of conservation 
of angular momentum of the particle, and it is not immediately obvious why it 
is inapplicable for short \ D . It is shown that there is a critical trajectory which 
is such that all trajectories inside it correspond to particles which are collected, 
whilst all trajectories outside it are those of particles which miss the probe. 
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Eq. (5) is then found to be valid if at the point of closest approach the critical 
trajectory grazes the probe surface, and ceases to be valid when the sheath field 
is so strong that the trajectory crosses the probe surface with a radial velocity 
component. This occurs when the distribution is such that the potential V at 
radius r satisfies the inequality 



p 


(?) 


Since in the absence of space charge the potential falls off less rapidly than (7), 
it follows that (5) is precisely applicable until some minimum value of \ D is 
reached which, unfortunately, has not been calculated though measured charac- 
teristics suggest that (5) is in fact applicable in the ionosphere for small values 
of probe radius. 

Thus we see that in principle the same information, i.e. the number density 
and mean energy of the particles, may be available from the accelerating as from 
the retarding characteristic, but that the form of the characteristic now depends 
both on the geometry of the probe and on \ D . Evidently such factors as the need 
for some means of supporting the probe will now result in the occurrence of er- 
rors which would be absent in the retarding regime. The advantages of the latter 
have caused it to be generally used. 

(c) Determination of Ion Temperature 

The method of determining T & will be evident from (a) above, but that used 
for T requires some discussion. By means of Eq. (4), the ion energy distribu- 
tion function can be obtained [Bowen et al (1964(a)), Boyd and Raitt (1964)] . 
Physically it is evident that the mean energy of each ion species of mass M + is 


10 




approximately E q = — M + v s 2 where v s is the spacecraft velocity, the> energy 
spectrum thus containing peaks corresponding to each of the ion species present. 
Then each peak will be broadened by the ion thermal motions. Medieus (1962) 
has shown that the energy distribution function for a single species is 



1 /^T /-(E-S 0 ) 2 -\ 

^ 2 ’ ttE o ex P ^4E kT L j 


( 9 ) 


when E q >> kT x , T l being the ion temperature. 

(9) shows that for the heavier ion species such as O* , the peak shape is nearly 
Gaussian with a width from which the temperature can be determined (Fig- 
ures 3 and 4). The mass separation of the major ions in, the F-region (H + , He + 
and 0 + ) is so large and the number densities of other constituents so small that 
the temperature determination is simple in practice. 

It is not essential to use a spherical geometry if the form of the distribution 
function is assumed so that only the parameters n Q and T are to be determined 
from it. The current voltage characteristic for a plane probe facing in the di- 
rection of motion of the spacecraft has been calculated from the one-dimensional 
equivalent of (8) by Whipple (1959), 

en o v s 

J p = 2 1 + erf « + V _ ^E L e_X 

o 
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where 



A curve fitting technique is then used to determine the concentrations of the 
major ion species and the temperatures from a measured characteristic. The 
procedure, which is much more difficult than for the spherical probe, has been 
described by Moss and Hyman (1968) and by Patterson (1969). 

(d) Examples of the Method 

Many different probe geometries have been used by different workers (see 
Boyd (1968) for a description). Here we shall briefly describe two examples of 
such experiments which are of particular interest on account of the use we shall 
make of data which has been obtained by them and because they' have been inter - 
compared and also compared with incoherent scatter results. 

The first is the cylindrical probe used by Spencer, Brace and their collabo- 
rators, which has been described for the Explorer XXXI satellite by Findlay and 
Brace (1969). The probe is a 0-057 cm. diameter wire 46 cm. long, with a guard 
cylinder over the wire at the support end to reduce the effect of irregularities in 
the electric field near the support. Such irregularities should not affect the 
evaluation of T e but may shorten the exponential region of the characteristic 
and lead to errors in the density. A sweep potential from -SV to +10.5V was 
applied to the probe with a repetition rate of 0.5 cycles/sec. The probe current 
was measured by two linear amplifiers each with two switched ranges of sensi- 
tivity, the outputs being telemetered. A generally similar arrangement has been 
used in rocket experiments except that the high telemetry sampling rate enables 


12 



a higher sweep rate to be used [Spencer, Brace, Carignan, Taeusch and Neumann 
(1965)] . Figure 1(c) shows the accelerating portion of the probe characteristic 
which is only roughly of the form given by (5) with a-= Since the current 
rises more rapidly than would be expected the difference from the expected form 
cannot be the result of sheath formation. It probably is due to the electric fields 
near the end of the probe remote from the guard electrode. Figure 1(b) shows 
the electron retardation characteristic, from which the method of determining 
the vehicle potential and also the probe current at space potential can be seen. 

The use of probes based on Druyvesteyn's relation (4) has been described 
by Bowen et al [1964(a), 1964(b)] and by Wrenn (1969). Small amplitude sine 
waves of two different frequencies in the audio range are applied to the probe, 

* I s 

and the intermodulation component, which is proportional to 3 2 j /3 V p 2 , is se- 
lected by a tuned amplifier. By applying a sweep voltage as well, the quantity 
_ i 

V 2 f (V) is traced out, from which if desired f (V) can be obtained. Figure 2 
shows a spherical ion probe. The central ion collector is 9 cm. in diameter. It 
is surrounded by an electroformed grid 10 cm. in diameter with holes 0.7 cm. 
diameter, the grid material being 65 /z thick. The grid transparency is 35%. The 
probe is mounted on a short boom to avoid the major effects of the spacecraft 
wake [see Section 3 (d)] . Figure 3 shows the probe characteristic for an ion 
mixture and Figure 4 shows a comparison between the observed and best-fitting 
characteristics . 

3. SOURCES OF ERROR IN PROBE MEASUREMENTS 

Probe measurements are potentially subject to errors from a number of 
sources. It will be shown in Section 6 that reliable results can be obtained 
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provided they are carefully taken into account. Here each of them will be ' 
briefly discussed. 

(a) Photo-emission 

If any part of the probe system is at a negative or small positive potential 

* 

with respect to the surrounding region, photo- electrons may be liberated by sun- 
light. Because of the large amount of the UV energy which is concentrated in the 
Lyman-a line at 1216 A, which has sufficient quantum energy to release electrons 
from all the materials likely to be used, the photo-current is not greatly depend- 
ent on the work function of the material chosen. 5 x 10 -9 amp/ cm" 2 is a typical 
value for tungsten, gold or rhodium. This is about equal to the positive ion com- 
ponent of the probe current for densities- in the range (0.3 - 3) x 10 5 cm" 3 so its 
effect must be eliminated in some way. Over the electron retardation character- 
istic the photo-current is constant so that it may be determined at a suitably 
large negative probe potential and simply subtracted at other negative potentials. 
Alternatively a suppressor grid can be used in front of the collector. Since the 
energies of the emitted photo-electrons do not exceed a few eV, a negative bias 
of this amount relative to the collector will suppress photo -emission from the 
collector although that from the grid still remains. When the probe is operated 
in the ac or Druyvesteyn mode, the effects of photo-emission are automatically 
suppressed over the electron retardation characteristic. The simplest solution 
of all, in some circumstances, is to ensure that the probe remains in the shadow 
of the spacecraft . 

(b) Geomagnetic Field 

In the theory briefly developed in Section 2, the effects of the geomagnetic 
field were neglected, which is entirely justified for the positive ions on account 
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of their large gyro-radius in comparison with the probe dimensions usually 
employed. For the electrons, where these quantities are of the same order, 
more careful consideration is needed. The presence of a magnetic field does 
not affect the electron energies, but by reducing ] o makes the condition (2) more ' 
stringent. It would therefore be expected that at sufficiently negative potentials, 

T e could still be measured by the use of (3), but that j o and hence n o could not 
be accurately measured as V p 0 . 

The application in the ionosphere of earlier theoretical work has been dis- 
cussed by Dote et al (1965). They show that the retarding potential method for 
T e remains valid at sufficiently great probe potentials, i.e. when 

V p > A B 2 d 2 (10) 

where B is the geomagnetic field intensity in oersted, d is the sheath thickness 
in cm. and the constant A is 8.8 x l(T 2 volts (oersted cm.) -2 for a plane probe 
and 2.2 x 10“ 2 volts (oersted cm.) -2 for a cylindrical probe. Since' in the iono- 
sphere B does not_ exceed 0.6 oersted and d is usually less than 2 cm, we find 
even for a plane probe that (10) requires V p > 0.126 volts which is not in prac- 
tice a very stringent condition except at low values of T g such as may occur in 
the D and E -regions. 

The requirement that the effect of the magnetic field be small at V p = 0 can 
be met only if the probe radius a is sufficiently small, in the case of a spherical 
probe less than 0.5 cm. in diameter. Smith (1969) has detected a magnetic field 
effect on a nose-tip probe maintained at an accelerating potential, the probe being 
approximately a cone with a base diameter of 2.5 cm. and an included angle of 11°. 
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Raitt (1969) on the other hand mounted two plane probes 2 cm. diameter with a 
guard ring of outer diameter 4 cm. on'ESRO IA. This spacecraft was stabilised 
with respect to the geomagnetic field and the probes were mounted so that the 
field was parallel to the surface of one' and perpendicular to that of the other . 

N 

The electron densities determined by the two probes agreed to within 10% except 
when one of the probes was in its own wake. These results are qualitatively at 
least in agreement with the conclusions of Dote et al. 

It appears that the geomagnetic field is unlikely to cause error in the 
measurement of T e in view of the practice of using an extended portion of the 
exponential characteristic in estimating its slope but for probes larger than 
1 cm', in dimensions, errors in the determination of n will occur. 

7 O 

(c) Collisions 

» 

V. 1 

The discussion of Section 2 will be invalidated if the charged particles being 
collected collide either with neutral or other charged particles near the probe. 

It is probable that the accelerating characteristic will be most affected. It is 
possible that a similar argument to that advanced in the case of the geomagnetic 
field can be used here also, and that at sufficiently negative potentials the meas- 
urement of T e is unaffected. The region of applicability of the Langmuir treat- 
ment in the ionosphere has never been precisely determined but it is usually 

J 

assumed to' apply when the mean free path exceeds either the sheath thickness 
or the effective probe dimensions, allowing for the effect of the electric field in 
the accelerating rdgime. This is at altitudes of 85 - 90 1cm for the dimensions 
usually employed. This assumption is of some importance4n view of an apparent 
discrepancy between incoherent scatter and probe measurements of T e below 
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130 km. which might be the result of unexpected collisional effects on the probe, 
though this does not seem very likely. 

Hoult (1965), Hoult and Ta-Jin Kuo (1966) and Sonin (1967) have treated the 
case of the probe where many collisions occur; it appears in this case that the 
charged particle number densities can still be obtained but information on the 
energy distribution is lost. 

(d) The Spacecraft Wake 

The formation of the wake has been the subject of many theoretical studies 
[for an account see Henderson and Samir (1967)] . Of these one which gives a 
simple physical picture and results which are in excellent agreement with obser- 

I 

vations is that due to Gurevich (1962). The relatively large energies of the ions 
in the spacecraft reference frame and the limitation of the spacecraft electric 
field to a fairly thin sheath together ensure that the ion trajectories are little 
affected by the spacecraft unless they are intercepted at its surface. The space- 
craft potential is usually attracting for the ions, and its velocity is comparable 
with or greater than the mean ion thermal speed even for rockets so that ions in 
the path of the spacecraft are simply swept up on the front surface, leaving behind 
a void. Since few ions have sufficient speed to overtake the spacecraft a conical 

y 

surface is formed as the ions behind the spacecraft move in towards the axis with 
the mean thermal velocity. Thus the wake consists predominantly of an approxi- 
mately conical rarefaction whose length is about a s v /C where a s is the 
spacecraft radius (assumed spherical). Since X. D < < a s , the electron density 
must now adjust itself to the ion density, except for a small excess required to 
charge the wake negative and so produce a density reduction in accordance 
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with (1). The wake thus acquires a negative potential of 5 - 10 kT e /e . The 
electron density in the wake, for a satellite at altitudes where the ions are pre- 
dominantly 0 + may fall to 1 % of the ambient value, the ion densities being even 
less. In front of the spacecraft, which reflects the bulk of the electrons, a 
smaller negative potential, less than kT e / e, will develop. This, acting on the 

ions, will slightly raise the density above the ambient value. Evidently, the space- 

/ 

craft significantly modifies the density distribution in its vicinity. Changes in T g , 
however, in the absence of significant collisions, are very much less likely to 
occur. Even the effects' of the density changes may largely be avoided by suitably 
positioning probes on a spacecraft which is stabilised relative to ihe velocity 
vector, or by the use of booms whose length is 1 - 2 spacecraft diameters. 

The Ariel I wake has been extensively investigated [ Samir and Willmore 
(1965), Henderson and Samir (1967)3 and shown to conform to this picture. 

Figure 5 shows the electron density depletion observed. 

(e) Spacecraft Potential 

The probe has so far been treated as an isolated electrode, but the potential 

distribution about it is actually influenced by the charge distribution on the space- 

* 

craft. Moreover, the spacecraft is used as the reference potential for the probe 
sweep circuits and the probe and spacecraft act as a double probe which as a 
whole can collect no net current from the ionosphere because its capacitance is 
small [Boyd (1968)]. 

\ 

The first of these effects will be greatest when the probe is mounted in the 
surface of the spacecraft. This case has been treated in detail by Whipple and 
Parker (1969) assuming that k D is very large so that space charge effects are 
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negligible. This is the worst case since the effect of space charge is to screen 

the probe from the influence of the spacecraft through the formation of a sheath. 

They show that, where the probe and spacecraft potentials have the same sign 

but the magnitude of the probe potential is smaller, or where the potentials are 

opposite in sign, a saddle point in potential develops in front of the probe. This 

\ 

is because the influence of the probe potential dominates near the probe, but that 
of the spacecraft at large distances. Thus for particles for which the spacecraft 
potential is repulsive — the electrons in the usual case — a potential barrier de- 
velops around the probe and the probe current density will not reach the value 
given by Eq. (3) if the total current drawn would exceed that across the barrier. 
In an extreme case the knee in the probe curve occurs at satellite potential in- 
stead of at space potential. This can be considered as the result of reducing j o 
m (2) and Langmuir recognised the possibility of such an effect in noting that a • 
probe placed in a re-entrant cavity might not be able to draw the full current. 
The probe characteristic for large potentials of the same sign as that on the 

i 

spacecraft is unaffected, so that in the usual case the electron retardation char- 
acteristic at large negative potentials is not affected, and measurements of T & 
can still be made. 

Results in agreement with these conclusions have been obtained by Wrenn 
(1969) on Explorer XXXI. The electron density measured by a plane probe 
mounted in the spacecraft skin is found to be much lower under some circum- 
stances than the ion density or the density obtained by the comparison topside 
sounder in Alonetto II. The difference is found (Figure 6) to decrease as X. D de- 
creases and the screening effect of the spacecraft sheath is felt. On the other 
hand T e measured by this probe is in agreement with that derived from other 
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probes on Explorer XXXI [see Section 6(a)] . The diameter of the guard ring of 
this probe was 4 cm. Another -plane probe 8.5 cm. in diameter mounted on this 
satellite [Donley (1969)] did not show this effect. 

The potential adopted by the spacecraft is" also of interest because this is in 
practice the reference for the probe sweep circuit. If for the moment it is sup- 
posed that the probe is at floating potential, i.e. the potential at which the.-net 
current is zero, then the spacecraft must also be at floating potential, since the 
capacitance is too small for a current to flow for any significant time [Boyd 
(1968)] . The electron current far exceeds that due to photo-emission or positive 
ions except at great heights so that the spacecraft must acquire a negative charge 

to reach floating potential, usually of about (5 - 10) kT e /e . Photo -emission will 

. -r 

tend to reduce or even to reverse this potential. Samir and Willmore (1966) have 
shown that the potential of Ariel I fits well this description. Whipple (1965) 
showed that certain spacecraft have taken up much larger potentials, invariably 
negative relative to the ionosphere. This was due to the presence of portions of 
the solar cell array and its wiring which were at large positive potentials rela- 
tive to the spacecraft, which, tending to draw a large electron current, drove the 
whole system negative. It was noted in Section 2(d) that space or spacecraft 
potential can be determined from the probe curve. Large negative values are 
only objectionable if the effect described by Donley and Parker occurs. The 
principal experimental problem, however, is that if the spacecraft attitude 
changes, as when for example it is rotating, then the different components of the. 
current may vary causing instability in spacecraft potential and distortion of the 
probe curve. This effect is most severe at low sweep rates and for ion probes 
which have an extended retarding characteristic. A.C. probes are almost immune 
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since only potential variations during the period of the lowest frequency used 
are important. 

Next consider what happens as the probe potential itself varies as a result 
of the applied ramp voltage. The spacecraft potential must also change so this 
must either be made sufficiently small or be accurately predictable if the change 
in probe potential relative to the plasma is to be known. No particular difficulty 
will be encountered with an ion probe as the current collected is always likely to 
be small, but the first condition can only be satisfied with an electron probe if 
its area is made 1000 - 2000 times less than the superficial area of the space- 
craft. If small parts of the spacecraft are at such positive potentials that they 

may approach space potential, their area should be subtracted from the per- 
^ * 

missible probe area calculated in this way. The alternative, rarely adopted, is 
to use a double probe in which the geometry and surface treatment of the space- 
craft are such that it, too, can he treated as a probe. This imposes rather severe 
constraints and the double probe system must usually be a separate ejected pack- 
age [Spencer et al (1962)] . 

(f) Variations in Work Function 

v. 

Also effectively in series with the probe sweep generator is any contact 
potential difference between the probe and the spacecraft, so that the stability of 
the work functions of each of these is of as much importance as that of spacecraft 

t 

potential. According to Wilson and Garside (1968), variations in work function 

♦ 

over the surface of a probe may vary from a few mV where considerable pre- 
cautions are taken to as much as 200 mV if the surface is handled. It is easy to 
show that when the probe is so negative that the entire surface is bn the 
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exponential portion of the characteristic /measurements of T e will be unaffected. 
The main eilect is to increase the founding of the characteristic near space po- 
tential and to introduce possible errors in the determination of n 0 . A potentially 
more serious problem [Sayers (1969)] are slow variations in probe work function 
during the sweep since these distort the probe characteristic. A.C. probes, due 
to the short effective measuring time, have an inherent advantage here also. 

(g) Aerodynamic Effects 

* 

It is usually assumed that aerodynamic effects will become important only 
in the D-region where the mean free path becomes smaller than typical rocket 
dimensions and a shock wave and its associated heating effects will develop. 

Smith et al (1968), however, noted a tendency for higher -electron temperatures 
tp be measured on the ascent than on the descent of the rocket trajectory even in 
the E-region. He used a nose-tip probe and presumably larger effects would be 
expected than with the. other geometries usually employed. In other measure- 
ments such a systematic assymetry does not seem evident. 

(h) Rocket Outgassmg 

If gas is emitted either from the interior of the spacecraft or from a spent 
or liquid fuel motor in the vicinity of the probe, the probe current may be 
reduced because of collisions between the charged and neutral particles and also 
because of the local depletion of ionisation as' the neutral gas cloud expands. 
These, effects can be controlled either by careful sealing or venting of the pay- 
load space and if necessary by ejecting the spent motor case. Since the evolution 
of gas will be greater on the ascent of a rocket trajectory the resulting errors 
will be greatest then. Comparison of ascent and descent results suggests that 
the effects on T e are not usually large. 
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4. MEASUREMENTS OF IONOSPHERIC PROPERTIES 
BY INCOHERENT SCATTER 

It was first suggested by Gordon in 1958 that a sensitive radar with a large 
aperture aerial might be used to study the ionosphere by irradiating it with a 
radio wave and observing as a function of height the intensity of the scattered 
wave resulting from Thomson scattering by the free electrons. Provided that the 
electrons are randomly distributed in space, the scattering cross-section per 
unit volume and solid angle is 

O' = n o ( n ) 

where cr c is the classical value for the electron cross section 



Moreover, on account of the thermal motions of the electrons Gordon believed 
that the scattered line would be Doppler broadened and that from its spectral 
distribution the electron temperature could be found. The first successful ob- 
servations were made at 41MHz by Bowles (1958) who found that the spectral 
width of the scattered line was approximately two orders less than this predic- 
tion. This had the important practical consequence that a substantially less 
sensitive system could be used to make successful measurements. As to the 
explanation of the narrow scattered line, Bowles suggested correctly that since 
the wavelength employed was considerably larger than fluctuations in the 

electron gas could not be considered without taking into account the effect of the 

' 

ions. He proposed that the effect of the Coulomb forces was to give the electron 
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density fluctuations (which are certainly predominantly responsible for the 
production of the scattered wave) the velocity spectrum of the ions, so that the 
spectral width was actually a function of T i . 

The problem of calculating the precise spectral distribution was tackled 
rather differently by Dougherty and Farley (1960), Fejer (1960), Salpeter 
[1960(a), 1960(b)] and Hagfors (1961) for the iso-thermal case, i.e. where T e = 

T , identical results being obtained. In each case, the spatial distribution of the 
electrons is considered random, so that the scattered power is proportional to 
the number density of electrons and the Lorentz method in which an electron is 
considered to oscillate sinusoidally in the electric field of the radio wave is 
used. This leads to the following expression for the scattering cross-section 
per unit volume of a small volume element V 

cr = cr c V < | n(K) | 2 > (13) 

where K is the difference of the incident and scattered wave vectors (K = An /X 
for backscattering), and n (K) is the spatial Fourier component of the electron 
density functions. < > represents as usual the expectation value in the sense of 
statistical mechanics. Physically this relation implies that the scattering results 
from refractive index variations of scale X/2. 

Where the methods differ is in the calculation of < [ n (K) | 2 >. Hagfors 

/ W\ , 

uses Boltzman's distribution function exp l — K^V ^ or a s ^te of the plasma of 
energy W. The plasma energy is the sum of that arising from the kinetic energy 
of the particles, and the electromagnetic energy. The latter term is attributed 
to electrostatic energy only, arising from the Coulomb interactions of the 
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particles, which can be expressed in the terms of the ion and electron density 
fluctuations giving 


< I n (K) | 2 > 

L 


n 1 + 2X 2 

° . p_ 

v 1 + 4X 2 

p 


(14) 


where 

p n U v 

Thus, as K - 0 (\ >> X D , the case usually encountered) cr -» n o a c /2. The 
scattering cross-section is one-half the value for Thomson scattering, the re- 
duction resulting from the effect of the ions in inhibiting fluctuations in n e . How- 
ever, as K -* °°, the scale of the fluctuations becomes small and the (coupling 
effect of the Coulomb fields negligible, so cr -> n Q a c as in (11). 

In order to obtain the spectrum of the scattered wave, the Fourier transform 
of the density fluctuations both in space and time, n (K, co) must be obtained. The 
component of frequency co will lead to a Doppler shift of co in the scattered wave. 
If the ion-electron mean free path were short compared with \ , fluctuations in 
density of wave vector K would be sinusoidal corresponding to an acoustic wave 
at a single frequency co = C ac • K where C ac is the acoustic velocity. When the 
mean free path is long (as in our case) such waves are heavily damped. Then, 
if the ions were absent, thes^e sinusoidal fluctuations would be replaced by ran- 
dom fluctuations resulting from the thermal motions of the electrons and it is 
clear that the apparent velocity of regions of high or low density will be related 
to the mean thermal speed of the electrons. Similarly, the spectrum of the ion 
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density fluctuations of wave vector K will be governed by the mean ion thermal 
speed. When Coulomb forces are considered the fluctuations in electron number 
density when KV D is small axe constrained to be nearly equal to those of the ions 
and take on correspondingly a frequency spectrum related to the ion thermal 
speed c . 

One aspect of this process is well revealed by the approach taken by 
Dougherty and his collaborators. They note that the density fluctuations (on the 
square of which the plasma energy depends in part) can be treated as a gener- 
alised coordinate of the plasma. Then, if a plasma admittance is defined in terms 
of the response of this co-ordinate to a sinusoidal driving force at the frequency 
co, the density fluctuations can be obtained from the real part of this admittance 
using the "generalised Nyquist theorem". If the admittance is not purely imagi- 
nary, then dissipation of the wave in the plasma must occur (corresponding of 

course to the generation of the scattered wave) and the question of the mechanism 

/ 

of this energy dissipation then, arises. In the absence of the wave, the electron 

density fluctuations in space will be random. When the radio wave is present, a 

degree of order is impressed on the electron gas which will now tend to fluctuate 

in density with wavelength X / 2 . Consider two regions of high density which are 

at some instant k/2 apart. Because of the thermal motions these will tend to 

move relative to each other but the influence of the wave is to maintain the sepa- 

/ 

ration at X/2. To do this, the electron velocities must be changed, i.e. work 
must be done on the plasma. This process is Landau damping and we see that 
the interaction of the waves is with the electrons whose kinetic velocity equals 
the "phase velocity" co/k of the fluctuations. 
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Figure 7 shows a normalised spectrum of the scattered line calculated by 


Hagfors (1961) for an atomic oxygen plasma. The variable 


1 



(15) 


is the ratio of the phase velocity to the most probable electron velocity. X p is 

the value of X when co = &> p , the plasma frequency. For X p << 1, the scattered 

line is a single peak of width X ~ 1 as expected from Thomson scattering. When 

i 

X p >> 1, the principal peak is reduced in width by a factor (M/m) 2 (M the ion 
mass, m the electron mass) which is 172 in this case. Some effect of the elec- 
tron motions can still be seen when X p is not too large, in the form of a weak 
extension of the spectrum to X = l.o. Note that, since the presence of the ions 
has narrowed the line considerably whilst reducing the total scattering only by a 
factor of two, the amplitude is much enhanced. A striking feature also to be ex- 
plained is the occurrence of a narrow scattered line, the "plasma" line, at X = X . 

Whilst acoustic waves are strongly damped in a plasma of low collision 
frequency, one collective motion still exhibited by the electrons is longitudinal 
plasma waves. In the absence of a magnetic field the dispersion relation for 
these waves is 


OJ - co p M + 3(477 X. D /X) : 
showing that for waves of wave vector K 


= % <1 + 3K 2 X d 2 ) 2 


a 




o> 

p 




(16) 


( 17 ) 
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The sinusoidal density fluctuations in these waves leads to a narrow line nearly 
at the plasma frequency. For a plasma in thermal equilibrium the power scat- 
tered in the plasma line is, relative to that in the main peak, K 2 which is much 
less than 1 [Perkins, Salpeter and Yngvessen (1965)] . However, the waves are 
as we have seen excited by electrons whose kinetic velocity is close to the phase 
velocity. For typical ionospheric values of co p and at wave frequencies of 
400 MHz, electrons of about 10 eV are required. Thus, the amplitude of the 
plasma line depends on the high energy tail of the electron energy distribution. 

We shall see in Section 7 that this is strongly enhanced in the daytime by a photo- 
electron flux resulting from photo-ionisation by sunlight. This enhancement is 
of the order of a factor of 50 and has been observed at Arecibo. The velocity 
distribution of the photo -electrons is strongly anisotropic as they are constrained 
to spiral up or down the geomagnetic field lines, so that the intensities of the two 
plasma lines on either side of the incident wave frequency are not necessarily 
equal, but enable the flux of photo-electrons towards and away from the trans- 
mitter to be determined. By observing the amplitude of the plasma line at vari- 
ous heights above 300 km. where the photo-electrons are free from collisions, 
the variation of <*> p enables the photo -electron spectrum to be explored from 
6 - 20 eV [Yngvessen and Perkins (1968)] . 

The above treatment has been extended to the case where T and T are 

6 1 

unequal, as is generally true in the ionosphere, by Fejer (1961) and Farley 
(1966). Results calculated by Loewenthal from Fejer's theory are shown in 
Figure 8. The spectrum now depends in detailed shape on T e /T , though the 
width is still controlled by 'T . 
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The influence of the geomagnetic field has been discussed by Farley et al 
(1961), Fejer (1966), Hagfors (1961). It is found to be surprisingly small until 
the propagation is very nearly (i.e. within about 5°) perpendicular to the field ' 
line. In the case of quasi-longitudinal propagation, by far the commonest case, 

i 

the relations already given in the absence of B can be accurately used. When 

the propagation is transverse however lines develop in the spectrum as shown 

\ 

in Figure 9 at multiples of the ion gyro-frequency. It was at first hoped that 
these would enable incoherent scatter to be used as a "magnetic mass spectrom- 
eter",. but the stringent observational criteria have prevented this hope from being 
realised. Physically, this is because the random motions of the ions along the 

field line (which are not inhibited by the field) smooth out the effects of their ro- 

% 

tation about the field except at very large angles, where the Doppler shift due to 
the movement along the field line becomes very small or zero. Hence themather 
narrow peaks due to the rotation about the field become smoothed out in the quasi- 
longitudinal case. 

For the special conditions of thermal equilibrium (i.e. where T e = TJ the 
fact that the Boltzman distribution is still applicable in a magnetic field (since 
the charged particle energies are not changed) enables the argument which led to 

i 

(14) still to be used. Thus, although the spectral distribution is modified by the 
magnetic field, the total scattering cross-section is unchanged. This is not true 
when T g f T i [Barron and Petriceks (1967)]. 

At altitudes near and below 120 km. the effects of collisions with the neutral 
particles must be considered [Dougherty and Farley (1963), Farley (1966),, 
Seasholtz and Tanenbaum (1969)] . Again if T e = T i , the scattered power is 
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unaffected because the Boltzman distribution may be used, but otherwise the 
scattered power, and in any case the spectrum, are changed.* 

The theory thus far described predicts that the scattered line will be sym- 
metrical about the incident wave frequency. If the ionosphere as a whole is 
moving relative to the transmitter along the line of sight, this will introduce a 
displacement of the scattered line. In this way the wind velocity in the line of 
sight can be measured. Such movements have been observed by Vasseur and 
Waldteufel (1969) who attributed them to the neutral atmosphere motion resulting 
from internal gravity waves. 

In addition to measurements of the scattered power and its spectral distri- 
bution,’ we can add that of the polarisation of the scattered wave, at least provided 
a rather low frequency (50 MHz) is used. The wave is subject to the usual 
Faraday rotation characteristic of quasi-longitudinal propagation in a magneto- 
ionic medium (Ratcliff e, 1956) and from the variation of rotation with height the 
electron density can be obtained. Many parameters — T g , T x , ionic composition 
[ Moor croft (1964)] , magnetic field (which can often be neglected), electron den- 
sity and collision frequency — influence the spectral shape and in the presence 
of experimental errors it is not possible to determine them all by curve fitting 
and to obtain a reasonable precision. Hence some simplifying assumption must 
usually be made, or one of these parameters must be independently determined. 
The scattered power is determined by T e /T i , electron density and the magnetic 
field (whose effect, once again, can usually be neglected). As it is difficult to 
determine the absolute scattered power with satisfactory precision, the usual 

*There is an interesting parallel between the results of this and the preceding paragraph and the 
corresponding conclusions in the case of the regarding regime of the Langmuir probe. 
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practice is to determine the shape of the electron density profile from the rela- 
tive variation of the scattered power and to normalise the profile at some point 
such as the F-region maximum by an independent observation of electron density. 
Alternatively, the process may be inverted if the electron density profile is known 
either from plasma line or Faraday rotation measurements [Farley (1967)] or 
otherwise [Maynard and Ducharme (1965)] and T e /T i determined. This de- 
termines one parameter which must otherwise be determined from the spectral 
distribution. 

The observed spectrum of the scattered line is also influenced by instrumental 
effects, since usually the operation as a radar requires that a pulse be transmitted 
in order to determine the range of the scattering region. Then, the scattered 
spectrum is the convolution of the pulse spectrum with that due to the ionosphere; 
this spectrum will be modified then due to the finite bandwidth of the receiver. 
These effects may be computed and the modified spectrum compared with that 
observed [Evans and Loewenthal (1964)] but rather more elegant procedures 
which involve either the use of a CW bistatic radar as at St. Santin - Nan§ay 
[Carru et al (1967)] or a double pulse techniques in which the auto-correlation 
function of the received pulse is determined [Farley et al (1967)] enable them to 
be avoided. With a pulsed radar, the resulting spectrum is an average over a 
height interval, typically 75 km., governed by the pulse length. With a CW sys- 
tem, the determining factor is the intersection of the two aerial radiation patterns 
and the height resolution may be as good as 5 km. 

Table I gives brief characteristics of some incoherent scatter stations. The 
reduction procedures differ somewhat at these stations, and so also do the num- 
ber of parameters determined from the scattered spectrum. For our purpose 
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Table I 

Characteristics of Some Incoherent Scatter Observatories 


Station 

Arecibo 

Jicamarca 

Millstone Hill 

St. Santin 

Latitude 

18 °N 

11.95°S 

42.6°N 

44.7°N 

Longitude 

67°W 

76.87°W 

71.5 °W 

2.2°E 

Radar type 

Pulse 

Pulse 

Pulse 

CW 

Aerial type 

Parabola 

Dipole 

array 

Parabola 


Aerial aperture, 
m 2 

7.1 x 10 4 

8.4 x 10 4 

1.9 x 10 3 

2 x 10 3 * 

Operating frequency 
MHz 

430 

49.92 

440 


Peak Transmitter 
Power MW 

2 

5 

t 


0.0751 

Pulse Length, /u. S 


200, 267 

500, 1000 

— 


* Receiving aerial; for the transmitting aerial the aperture is 8.10 3 m 2 
tMeaa transmitter power. 


it is sufficient to note that T e and T i are determined separately at each, except 
that the assumption may be made m the E-region that T e = T t . 

5. INDIRECT DETERMINATIONS OF PLASMA TEMPERATURE 

Many attempts have been made to determine the plasma temperature by 
methods less direct than those described in the previous sections, e.g. from the 
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altitude variation of ion composition. The method most extensively employed 
uses observations of the altitude profile of electron number density in the topside 
ionosphere. The scale height H defined by 

B(ln n o )/3Z = - i/H (18) 

is assumed to have the value appropriate to an isothermal ionosphere 

H = 2kT p /M t g (19) 

where k is Boltzmann's constant, g the local value of the acceleration of gravity, 
the mean positive ion mass and T p = ( T e + T x ) / 2 the plasma temperature. 
Either H is determined for an altitude at which it is believed that some con- 
stituent such as O* is so predominant as to give M i as the mass of that con- 
stituent, or an attempt is made to determine an ion composition which will re- 
produce the observed number density profile. 

The equilibrium of the ionisation in the topside F-region has been discussed 
by Bauer (1969). The equation of continuity 

c>n o /Bt = q - L + V • (n o v) (20) 

where, v is the transport velocity of the ionisation and q and L are respectively 
the volume production and loss rates of ionisation, shows that several scale 
heights above the F-region maximum where q and L have become small, and 
under steady state conditions where Bn Q /3 t approaches zero, 

V • (n 0 v) = 0 (21) 

This has two solutions, n o v = constant and v = 0, the electron density pro- 
files for which will approach one another at increasing altitudes as the 
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ambipolar diffusion coefficient increases. The solution v = 0 leads to 

2T 

H = 2 3T p /3Z + M ig /k (22) 

In the presence of the geomagnetic field this result will still be correct provided 
that it is applied to the electron number density variation with altitude measured 
along the field line since this is the direction in which the ionisation diffuses. 

The importance of this is greatest at low dip latitudes [Rishbeth et al (1966)] . 

The first point which arises is the applicability of (22) in view of the as- 
sumptions, particularly that v = o, which have been made in its derivation. For 
this simultaneous measurements of n o , T p and ion composition are required. 
Brinton et al (1969) and Brace et al (1969) analysed profiles obtained in a rocket 
flight and found the scale height between 300 and 400 km to be only about half the 
value given by (22), increasing between 400 and 500 km to the expected value 
without any large corresponding change in T . The difference was attributed to 
ionisation movements, the combined results of an upward flux of H + and a neutral 
wind. Wrenn and Smith (1969) obtained scale height values from topside ’sounder 
electron density profiles from Alouette II and values of M x and T p from a.c, 
probes on Explorer XXXE for effectively identical locations and deduced the re- 
quired values of 3T p /3Z from (22). Values ranging from -1 to +1.5 K/km were 
obtained. Colin et al (1969), in a very similar analysis using different data from 
the same spacecraft, obtained similar results. Gradients of 0.2-1. 0 °K/km are 
also obtained from incoherent scatter profiles of T e and T at lower altitudes, 
and values of 0.6-1. 2 °K/km maybe deduced from satellite observations of T e at 
different altitudes [Prasad (1968), Moor croft (1969), Willmore (1965), Mahajan 


34 



and Brace (1969) J , except that negative gradients have not been observed. Over 
a limited altitude range negative gradients in T i could result from composition 
changes [Banks (1966b)] resulting in negative gradients of T . It seems likely 
that the temperature gradients deduced from the diffusive equilibrium density 
distribution thus agree well with those measured directly, hence confirming the 
applicability of the scale height expression over a wide range of ionospheric 
conditions above 500 km. at least to an accuracy of ±10%. 

j 

On the other hand, the application of this result to the deduction of T p is 
very difficult, since neither the temperature gradient nor the proportion of light 
ions is easily predicted. Matuura and Ondoh (1967) analysed Alouette n profiles 
and showed that, even if the atmosphere was assumed isothermal, it was impos- 
sible to determine whether ox not a small proportion of He + was present, leading 
to uncertainty of about 20% in T p from this cause alone. Taieb (1969) combined 
Alouette electron density profiles with near-simultaneous values of T p obtained 
from incoherent scatter, and so was able to deduce values of the temperature 
gradient and the light ion proportion, but it seems unlikely that this would have 
been possible unless T p had been known. In all other attempts to deduce T pJ the 
temperature gradient and the effects of contamination by light ions were simply 
ignored. This, for results obtained near 500 km, may lead to errors of between 
-25% and +100%. Evidently these results are less reliable than those obtained 
by the other methods we have described, and we shall not therefore discuss them 
further. 
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6. COMPARISON OF TEMPERATURES OBTAINED BY 
PROBES AND INCOHERENT SCATTER 

/ 

Measurements of temperature obtained by probes mounted on satellites and 
by incoherent scatter .stations are somewhat complimentary, the satellites giving 
good geographical coverage and also results to higher altitudes, whilst the radars 
give good seasonal, diurnal and altitude coverage up to their maximum altitude 
of operation. It is therefore of importance to know how far both sets of obser- 
vations are consistent and how well different variants of each basic method can 
be compared, though the latter comparison is possible only for the probes. The 
first attempt at such a comparison was made by Evans (1965c) using incoherent 
scatter results from Millstone Hill and probe results from Explorer XVII and 
Ariel I. Large day-by-day fluctuations in T e occur at any location, so that un- 
less strictly simultaneous observations for a given location are used, a suf- 
ficiently large number of observations to give accurate mean values must be 
used. Moreover, the sets of observations used by Evans did not cover the same 
epoch. Evans did make one near-simultaneous set of observations with a rocket 
launch of the cylindrical probe, but rather large scatter in Ms results makes the 
comparison inconclusive. Subsequently several sets of observations have been 
obtained which enable a more satisfactory comparison to be made. 

(a) Intercomparison of Langmuir Probes 

In the ISIS series of satellites, Explorer XXXI was launched as a piggy-back 
to Alouette n. The. latter carried a topside sounder and a cylindrical probe, 
whilst the former carried cylindrical and planar Langmuir probes operated in 
a dc mode, and a plane probe operated in an ac mode [Brace and Findlay (1969), 
Donley (1969), Wrenn (1969)] together with a plane ion trap operated in a dc mode 
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and a spherical ion probe operated in an ac mode [Donley (1969), Wrenn (1969)] . 
The cylindrical probes on the two spacecraft were similar and so afford an oppor- 
tunity of assessing the possible influence of the spacecraft on the probe. The 
two spacecraft separated at 13 km/day, passing through the same region of 
space with a time delay which increased from zero at launch at a rate of 1-6 sec/ 
day. Brace and Findlay (1969) compared corresponding probe curves and con- 
cluded that they were generally in agreement to better than 10%, nor was there 
evidence for any systematic difference between the results for the two space- 
craft, despite the presence. on Alouette II of long aerials and the topside sounder 
transmitter. 

The three Langmuir probes on Explorer XXXI were inter compared on five 
passes where T g ranged from 1500-4000°K [Donley et al (1969)] . The results 
were generally consistent to within 20%. The ac mode Langmuir probe gave a 
result which on'an average was 12% less than the cylindrical probe, whilst the 
dc mode plane probe was 3% lower than the cylinder. It appears that there is no 
evidence for an effect of geometry on the result, but the ac technique gives a 
result 10-12% lower than the dc version. No intercomparison of the probe cali- 
brations was carried out prior to launch, so any effects of calibration errors are 
also reflected in these results. The rms deviations*of-the ratios of the tempera- 
tures obtained by the two pairs of probes were 13% and 9% so that considerable 
non-systematic differences also existed. Since the probe measurements were 
not strictly simultaneous, these could arise from temperature fluctuations along 
the satellite track, but examination of the results shows that the temperature 

ratios are usually more consistent within each pass than they are from pass to 

v 

pass, which suggests that this is not the explanation. The discrepancies are 
greatest at the lowest altitudes. 
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The ion temperatures cover a range of 1000-3 000°K. The results obtained 
on the five passes, on three of which H + was the dominant ion and on two 0 + , 
show marked discrepancies, the mean ratio of T x for the plane to the spherical 
probe being 1.454 with no clear dependence on ion species. The latter feature is 
surprising, since the most likely source of. error seems to be stray electric 
fields in the sheath or near the grids used, which should affect H + far more than 
0 + on account of the much lower mean energy. Knudsen (1966) concluded however 
that these errors would actually be negligible for a plane probe and it is not obvi- 
ous that they should be greater for a spherical probe. The spherical probe is 
operated in an ac mode which we have seen also gives lower electron tempera- 
tures. Donley et al (1969) concluded that the plane probe T values are in better 
general agreement with incoherent.scatter results which were however not taken 
simultaneously. Subsequently simultaneous results were obtained for the St. Santin 
and Malvern incoherent scatter stations which provide strong support for the 
spherical probe results. These will be described in the next section. The unex- 
plained discrepancy between the two probes is clearly unsatisfactory and suggests 
that measurements of T x were not made by these methods with the same relia- 
bility as T e . 

A comparison of the values of n fi with those obtained from the Alouette II 
topside sounder was also made. If the results from the spherical ion probe 
which was mounted on a boom are ignored, a small systematic discrepancy re- 
mains, the probe values of n e being about 3% higher than those given by the 
topside sounder. Since the probe measurements are all made on the forward 
side of the satellite relative to the velocity vector, this may be the density en- 
hancement there due to the ram effect. 
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(b) Comparison of Probes with Incoherent Scatter 

Tlie first attempt to compare probe and incoherent scatter observations was 
that of Hanson et al (1969) using passes of Explorer XXXII near the incoherent 
scatter radar at Jicamarca. A comparison of 5 passes distributed in local time 
showed the probe results for T e to be systematically higher, on an average by a 
factor of 1.7. On 1968 March 17, rockets carrying (inter alia) cylindrical probes 
were launched near the Arecibo station in Puerto Kico at local times of 03 h 30 m 

, v 

and 15 h 30 ra [Brace, Carlson and Mahajan (1969)] . The day profiles of T e show 
agreement within 5%, whilst at night the probe results were systematically 10-15% 
the higher, though the successive radar profiles show fluctuations of this order 
which may account for the difference. The altitude range of comparison was 
130-290 km in the day and 210-290 km at night. 

. Booker and Smith (1969) have compared cylindrical probe results (T s ) from 
Explorers XXXI and XXXII and ALouette II with incoherent scatter values (T B ) of 
T e taken near -simultaneously at Jicamarca, Arecibo and Millstone Hill and have 
concluded that the relation between these is 

T b = 0 84 T s - 400 (23) 

This corresponds to a considerably 'larger disagreement between the two esti- 
mates of T e than in the case of the rocket flights just mentioned. 

On 1968 April 23 a day pass of Explorer XXXI near St. Santin enabled a' 
comparison of the ac probe and incoherent scatter results to be obtained 
[Henderson et al (1969)] . Unfortunately the radar profiles could not be taken 
above 425 km whilst the satellite altitude was 575 km at closest approach. 
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Figure 11 shows the results from which it is clear that the satellite observations 
represent reasonable extrapolations of the radar data with a gradient of 0.8°K/km 
for and 2.0°K/km for T e . A discrepancy exceeding 200°K or 10% for T e and 
15% for would be unlikely. In particular, raising by a factor of 1.454, as 
required by the ion trap [see Section 6{a)] would bring T l = T e at 600 km, which 
also seems unlikely. 

Similar (but more extensive) comparisons between the Explorer XXXI ac 

probes and Malvern incoherent scatter values of T e and T i were made by Taylor 

and Wrenn (1970). They found the values of T e over the range 1000-3000°K to be 

in excellent agreement. The values of T i over the range 700-1500°K showed a 

large scatter but there was little evidence of systematic disagreement. 

* 

The difference between the ac and cylindrical probes on Explorer XXXI may 
either reflect the characteristics of the particular instruments on that spacecraft, 
or a more fundamental effect of the different techniques used. This point is of 

r 

sufficient interest to' suggest a more extensive comparison of the earlier satellite 
measurements of Ariel I and Explorer XVII with the radar data. No direct com- 
parisons were of course made. The principal difficulty arises from the day-to- 
day variability of electron temperature. This will be minimized in observations 
made at moderate altitudes near the geomagnetic equator since here T e never 
rises far above T n . Figure 11(a) shows that Jieamarca and Ariel I values of 
agree well with determinations of T n from the drag on Explorer I, and fairly well 
with exosphere temperatures estimated from the Harris and Priester formula 
[Priester, Roemer and Volland (1967), Jacehia and Slowey (1966)] . The Ex- 
plorer XVII T e values, however, agree well only after adjustment by Eq. (23). 
These observations were made at a geomagnetic latitude of 10°N,''So that in the 
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day temperatures up to 250°K above the equatorial values might be expected. The 
adjusted values are 100~150°K higher. 

A similar comparison can be made for latitudes near 40°N geomagnetic 
[Figure ll<b)] . At night the Ariel I and adjusted Explorer XVII observations are 
in good agreement with Arecibo values. The pronounced seasonal variation which 
is a feature of the Millstone Hill temperatures [see Section 7(c)] makes com- 
parison with that station more difficult, but the summer temperatures do not 

agree too badly. It seems possible that the Ariel I temperatures are about 100°K 

\ 

low. During the day, the variability of T g — indicated by the range shown for the 

Arecibo values — makes the comparison no longer possible. 

1 

From these results it seems most probable that the Ariel I ac probe gave 
results close to the radar values and hence that the discrepancy between the two 
types of probes, reflects a genuine difference resulting from the different tech- 
niques employed. 

i 

We can conclude that the cylindrical and planar dc proves are in good agree- 
ment, that the de probes give higher values of T e than either the ac probes or 
the incoherent scatter, and that the last two techniques are in good agreement. 

The high values obtained by the dc probes are not associated with the geometry, 
the use of retarding grids or the Spacecraft, and presumably therefore arise as 
a result of measuring the dc current rather than the derivatives of the charac- 
teristic. The high consistency obtained between the cylindrical probes on different 
spacecraft shows that this is not connected witlfthe current measuring system it- 
self. The most likely explanation appears to be that advanced by Sayers (1969) that 
the effective contact potential of the probe surface varies slowly with the probe 
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potential, the variations being small in the rather short time period of the low 
frequency waves used in the ac technique (2mSec) but significant in the rather 
longer time (tens to hundreds of mSec) required to sweep through the retarding 
characteristic. Further laboratory and ionospheric experiments would be needed 
to confirm this hypothesis. 

In the E-region, , below 130 km, there is no such direct comparison data 
available. The situation there will be further discussed in Section 8. 

7. TEMPERATURE OF THE MIDLATITUDE IONOSPHERE ABOVE 140 km 

It will be convenient to consider aspects of the behaviour of the temperature 
distribution in the ionosphere under several different headings. Above 140 km 
the major and probably the only source of energy heating the ionosphere is solar 
radiation, except under conditions of exceptional disturbance or at high latitude. 
We can thus divide the ionosphere into regions above and below 140 km. Under 
very disturbed conditions, such as in the vicinity of aurorae, this division ceases 
to have any significance and similar processes occur at all altitudes. Below 
140 km, cooling of the ionosphere by collisions becomes very efficient, and it is 
very difficult to account for the temperature enhancement often observed. 

* 

The first account of the processes determining the ion and electron tempera- 
tures in the region with which we are now concerned was given by Hanson and 
Johnson (1961) and was soon extended by Hanson (1962). When a neutral particle 
is ionised by a solar EUV or X-ray photon, the photo-electron^produced has as 
kinetic energy the excess of the quantum over the ionisation energy of the neutral. 
This energy, which is typically 10 eV , is then imparted by collisions to the 
neutral gas and to the ambient electrons, the latter predominating at high altitudes 
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and when the photo-electron energy falls below 1.5 eV. Thus, out of the initial 
energy, at 300 km, most of it is passed directly to the ambient electrons, the 
amount falling at lower altitude's to 1 eV at 100 km. This energy is rapidly 
randomised by electron-electron collisions, resulting in an enhancement of T e - 
above the gas temperature T . Because of this temperature difference, the 
electron gas now loses energy by collisions both to the positive ions and to the 
neutral gases. This results in an enhancement of T , so that the positive ions 
also lose energy to the neutral gas in ion neutral collisions, the ions taking up 
an intermediate temperature between T e and T g . 

Thus all the photo-electron energy terminates in the neutral gas and is 
finally removed by downward conduction and by radiation. The conductivity of 
the atmosphere is high so that through most of the region T g varies little with 

i 

altitude [Figure 12(a)] and its geographical variation is largely determined 
simply by the zenith distance of the sun [Figure 12(b)] leading-to a small de- 
crease in T g at high latitudes. These figures are taken from the review by 
Jacchia (1965) which discusses other influences on T g . 

Above about 300 km, the range of the photo-electrons emitted in an upward 
direction becomes very long, so that an upward flux through the magnetosphere 
and into the conjugate ionosphere is produced. Some energy is deposited in the 
magnetosphere, raising the temperature there, whilst the rest contributes to 
heatmg the conjugate ionosphere. The magnetospheric heat is removed by- 
thermal conduction to low altitudes where it is passed to the neutral gas by col- 
lisions so that a gradient of T e and T along the field lines is established 
throughout the magnetosphere and upper ionosphere. 
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This basic picture has undergone a good deal of refinement and elaboration 
without modification in essentials. In subsequent sections a detailed account will 
be given and compared with observation. Before* doing so, however, it will be 
natural to ask whether -under conditions which depart so markedly from thermal 
equilibrium, it is permissible to talk of kinetic temperature at all. The photo- 
electrons of course represent an enhanced high energy tail in the Maxwellian 
distribution which was shown by Hoegy et al (1965) to contain 10" 3 - 10" 4 of all 
the electrons present. The effect of an electric field component perpendicular 
to the geomagnetic field was considered by Megill and Cahn (1964) who showed 
that for a field of 100 mV /m, an extremely large value, a depletion of the tail 
occurs above 1.6 eV, the relative deficiency amounting to about 10" 3 of all the 
particles. Dalgarno and Henry (1965) considered the possible effect of photo - 
detachment from O a _ in the D-region which has a similar effect to ionisation at 
greater heights in releasing energetic electrons, and found that the enhanced tail 
included about 10" 3 of the electrons. Thus it appears likely that electron-electron 
collisions are effective at maintaining a Maxwellian distribution except above 
1 eV where deviations involving about 10" 3 of the electrons are to be expected. 
Banks (1968) has considered the effects of collisions on the ions. For an H + -0 + 
mixture, in which H + is in the minority, he finds that a Maxwellian distribution 
will be maintained for 0 + , but significant distortion of the velocity distribution 
for H + will occur. Thus when light ions are minor constituents we may expect 
them to have non- Maxwellian distributions, although the major ions will still 
have a Maxwellian distribution. Except in this case therefore, no difficulty is 
likely to arise in the interpretation of experimental results, or in the idealisation 
of the velocity distribution in theoretical computations by the use of a Maxwellian 
function. 
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Hitherto, the linearity of the semi-log plots obtained with Langmuir probes 
has been considered strong evidence for the existence of a Maxwellian energy 
distribution for the electrons . On the other hand, the discrepancies between the 
dc probes and incoherent scatter could result from non-Maxwellian distributions 
at low energies, to which the semi -log plots are not very sensitive. This possi- 
bility should be borne in mind in the subsequent discussion. 

(a) Energy Balance During The Day 

All solar radiation at shorter wavelengths than the photo-ionisation thresholds 
of the atmospheric gases contributes to the energy input in the day. Of the excess 
quantum energy of the ionising photon which appears as kinetic energy of the 
photo-electron, a fraction only, called the heating efficiency e , is available to 
heat the electron gas. The remainder is lost in collisions with other particles, 
predominantly those of the neutral gas. The starting point in calculations of e 
and the total heat input Q (which below the photo-electron escape level is e 
multiplied by the rate of ionisation) is a measurement of the absolute intensity 
of the solar spectrum [see, e.g. Hinteregger (1969)] . Together with a model 
atmosphere and laboratory or theoretical data on photo-ionisation cross-sections, 
this enables the rate of production of photo-electrons of variQus energies to be 
determined. Several collisional processes contribute to the degradation of this 
initial energy. At higher altitudes, where the neutral particle concentration is 
sufficiently low, Coulomb elastic collisions with the ambient electrons dominate, 
and nearly the entire photo-electron energy is available for heating. Here e 
approaches the mean initial kinetic energy. Lower down scattering by the neutral 
gas first predominates. For photo-electron energies above 5 eV , the important 
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colLisional process is electronic excitation-of O, 0 2 or N 2 , 'but for energies 
between 1.5 eV and 5 eV this is replaced by vibrational excitation of N 2 . Only 
when vibrational excitation becomes ineffective at still lower energies is the 
energy loss to the electrons once again important, so that only the last 1.5 eV of 
the photo-electron energy contributes to the heating and e = 1.5 at these altitudes, 
e and Q have been calculated in this way by several workers [Hanson (1962), 
Dalgarno et al (1963), Dalgarno et al (1967), Dalgarno et al (1968b), Herman and 
Chandra (1969), Nagy et al (1969)] , with typical results shown in Figure 12. In 
these calculations the escape of the photo-electrons is neglected so that it must 
be remembered that not all of the available energy is deposited where it is re- 

i 

leased above 300 km, although e correctly indicates the relative importance of 

the different scattering processes at these altitudes. Since both the production 

and loss rates of the photo-electrons are obtained, the number density distribution 

may also be calculated [Dalgarno et al (1963), Hoegy et al (1965), Nisbet (1968)] , 

(Figure 12). At lower altitudes the structure in the energy spectrum resulting 

from the transition from one collision process to another can be clearly seen; 

* 

this disappears at greater heights where only collisions with the ambient elec- 
trons are important. At these altitudes the spectrum is rather flat, extending 
from about 1 eV to 30 eV, The energy spectrum of the escaping photo-electron 
flux will be similar to this. ^ 

The values of Q shown in Figure 12 neglect the effect of photo -electron 
excape. From an analysis of Ariel I data which showed high values of T e to be 
maintained to 1200 Ian, Willmore (1964) noted the need for a source of energy at 
this height in excess of the direct solar input and suggested Coulomb collisions 
with the ambient electrons of the photo -electrons which were spiralling upwards 
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from lower altitudes as Hanson proposed. Geisler and Bowhill (1965b) called 

this non-local heating and they, from a discussion of the escape process, showed 

/ 

(Figure 15) that it exceeded the local heating by an order of magnitude a few 
scale heights above the escape level, which is at approximately 300 km. Physi- 
cally it is evident that the escape flux will be a fraction of all the ionisations 
occurring above the escape level z e , the fraction depending on the proportion of 
electrons emitted in an upward direction, and on the proportion of these which 
are scattered downwards. Geisler and Bowhill found for the escape flux G the 
expression 

HQ(z.) 

0 = -i— (24) 

where H is the scale height of the ionisable constituent and Q(z e ) the ionisation 
rate at the escape level. A detailed discussion of the process is complicated. 
The angular distribution of the emitted photo-electrons must first be taken into 
account as well as the effects of the magnetic field along which they spiral 
[Mariani (1964)] . Scattering by both neutral particles and electrons will be im- 
portant in defining the escape level [Geisler and Bowhill (1965b)] . A single 
collision will not always be sufficient to stop an electron moving upwards and 
multiple scattering must be taken into account [Geisler and Bowhill (1965b), 
Lejeune and Petit (1969)] . Nisbet (1968) has shown that in estimating the escape 
flux at various phases of the solar cycle, variations in the escape level due to 
changes in atmospheric temperature are important. He finds values for the 
particle and energy fluxes as shown in Table II. 

These results can be compared with several measurements of the photo- 
electron flux. It was noted in Section 4 that the upward and downward fluxes 
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Table II 


10.7 cm Solar Flux 
(xlO -22 W/cm 2 sec Hz) 

Escaping 

Photo-Electron Flux 
(el/cm 2 sec) 

Escaping 
Energy Flux 
(eV/cm 2 sec) 

75 

5.24 x 10 s 

5.43 x 10 9 

150 

4.81 x 10 8 

5.70 x 10 9 

225 

7.28 x 10 s 

9.81 x 10 9 


can be deduced from the intensity of the plasma line in incoherent scatter obser- 
vations. Moreover, a portion of the electron energy spectrum can be scanned 
from the height variation of the line intensity since as the height and hence the 
plasma frequency vary, so does the electron velocity required to excite the 
plasma waves. Observations of this kind made at Areeibo by Yngvesson and 
Perkins (1968) showed the spectrum to fall smoothly from 6 to 20 eV. In the 
day, upward and downward fluxes were about equal, showing that at this latitude 
most photo-electrons survived the passage through the magnetosphere. In the 
period 1966 Dec-1967 Mar. soon after the solar minimum, the particle flux was 
in the range (3-10) x 10 8 el/em 2 sec and the energy flux (3-13) x 10 9 eV/em 2 sec. 
Measurements made in 1965 Dec. -1966 June by Narasinga Rao and Donley (1969) 
and by Maier (1969) using retarding potential analysers on Explorer XXXI show 
the spectrum to have a broad peak from 3-13 eV, the maximum particle fluxes 
being 2.5 x 10 8 el/cm 2 sec. All these measurements are in good general agree- 
ment with the figures in Table II and with the calculations of Nagy et al (1969) 
but are substantially greater than those calculated earlier by Geisler and Bowhill. 
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Once the photo-electron flux and spectrum are known the energy gained by 
the ambient electrons is readily calculated from the expression given by Butler 
and Buckingham (1962) for the energy loss by Coulomb collisions by electrons 
of energy E, which reduces when E > 3 eV [Dalgarno et al (1963)] to 

dE 

-^ = - 1 95 x 10" 12 n e /E eV/cm 2 (25) 

A number of collisional processes contribute to the cooling of the electron 
gas by the neutral atmosphere, including elastic collisions, rotational and vibra- 
tional excitation of 0 2 and N 2 and excitation of 0( 3 P) to the *D level [Dalgarno et al 
(1967b)] . However, Dalgarno and Degges (1968) showed that transitions by elec- 
tron impact amongst the three fine structure levels of 0 

e + 0 ( 3 Pj) - e + 0 ( 3 P J1 ') 

where it was assumed that collisions with the neutral gas were sufficiently rapid 
to maintain the distribution amongst these levels close to that corresponding to 
the gas temperature T , would represent the dominant loss process to the neutral 
atmosphere. The relative rates of these processes are shown in Figure 16*. 

Above 300 km, the falling density of 0 makes this process less effective 
than the direct energy transfer to the positive ions in elastic collisions. This 
leads to a loss rate 

L ie = t^ 3/2 (T e " T i > eV/cC SeC < 26 > 

i e 

'In the discussions of this problem published prior to that of Dalgarno and Degges, and in some 
which appeared subsequently, the neglect of fine structure cooling leads to a serious underesti- 
mate of the electron cooling rate. 
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[Banks (1966c)] , where T x and n i are the temperature, and number density of 
ions of mass (in amu). Because of the more efficient energy transfer to light 
ions the thermal coupling of the ion and electron gases increases with altitude. 

The processes by which the energy thus gained by the ions from the electrons 
is communicated to the neutrals have been extensively discussed by Banks (1966c). 
The energy transfer rate is 

L . s = A ", CT - T,) (2Y) 

where n g is the number density of a neutral constituent. For ions which are not 
moving in their parent gas (e.g., He + in O) the energy loss is largely due to 
elastic collisions and may be calculated from the polarisability of the gas . The 
values of A p given in Table III have been obtained in this way. The errors re- 
sulting from the use of the polarisation interaction in the ease of H + in He at 
temperatures higher than 300°K have been considered by Dalgarno and Dickinson 
(1968) who whow that the error does not exceed 20%. For ions which are moving 
in their parent gas (e.g., 0 + in O) an important additional contribution to the 
energy loss arises from charge exchange which has a large cross-section at 
low energies in this case. It is also an effective cooling process as the initial 
velocities of ion and neutral particle are inter charged in the collision. The effect 
is important when (T^ + T ) exceeds 500-1000°K, i.e. generally in the atmosphere. 
The appropriate loss coefficient A cg is given also in Table HE. 

The steady state value of the ion temperature will now be established be- 

* 

tween T e and T g at such a value that the difference between the energies gained 
and lost by collisions with the electrons and the neutrals is ]ust balanced by the 
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Table III 


Energy Transfer Coefficients (eV/cc sec) For Ion-Neutral Collisions 

According To Banks (1966c) 


Particles 

A p x 10 14 

Particles 

A ce (T, + X 10 15 

o + -n 2 

6.6 

H + - H 

14.0 

0 
+ 

1 

o 

KJ 

5.8 

0 + - O 

2.1 

0 + - He 

2.8 

N + - N 

2.1 

0 + - H 

3.3 

He + - He 

4.0 

H + - O 

3.5 

CN 

O 

i 

+ <N 

o 

1.4 

H + - N 2 

3.1 

n 2 + -n 2 

2.7 

H + - He 

5.5 



He + - 0 

5.8 



He + - H 

10.0 



He + - N 2 

5.3 




net transport due to, thermal conduction. If for the moment the latter contribu- 
tion is ignored, this requires that L = L, . At low altitudes where n is large, 
~ T g ; at greater altitudes (in practice at 600-800 km) T x ~ T e . Because of 
the dependence of L ie on M i ~ 1 , and because of the different constants which enter 
into L. the different ion species present will tend to take up different tempera- 
tures [Banks (1967), Dalgarno and Walker (1967)] . The energy transfer rates 
between different ion species are about an order greater than those to the neu- 
trals, so the variation in T i is substantially less than T e - T . Moreover, the 
collision rate between ions of the same species is larger yet except for the 
minor constituents so that as we have already noted a Maxwellian velocity 
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distribution is maintained except for minor species. Banks finds that T(H + ) - 
T(0 + ) may reach 170°K at 310 km (or a rather larger value when the distortion 
of the velocity distribution of H + is considered) whilst the difference for He + is 
only one-quarter as large. 

When T ~ T g at considerable altitudes the effective energy loss rate for 
the electron gas becomes L ig , which falls in proportion to n e n g , whilst from 
(25) the energy input due to the photo-electrons is proportional to n . Thus at 
some altitude due to the fall in n g , we might expect very large values of T e . In 
fact the temperature rise in the topside ionosphere is limited by the fairly high 
conductivity K of the electron gas along the magnetic field lines [Hanson (1963), 
Willmore (1964), Geisler and Bowhill (1965b)] . Spitzer (1956) gives 

K = 20(2/tt) 3 / 2 (kT e ) 5/2 k/m 1 / 2 e 4 In A 

= 7.7 x 10 s T e 5 / 2 eV / cm sec deg K 

where A is a cut-off parameter for the Coulomb interaction, taken to be 15 in 
the ionosphere (may be 13-19). Above 500 km a temperature gradient of about 
l-2°K/km is established which is sufficient for all the energy deposited to be 
conducted down to altitudes where it can be passed to the neutral gas by colli- 
sions. Whereas in the lower ionosphere the temperature is governed by the 
balance between local heating and collisional cooling, in the upper ionosphere it 
is determined by non-local heating and thermal conduction. 

This situation was first discussed by Geisler and Bowhill (1965b) and sub- 
sequently by Nisbet (1968). It is found that some effect of thermal conduction is 
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observed even near and below the F-region peak in smoothing out local irregu- 
larities of temperature which might occur. Geisler and Bowhill also suggested 
that it may cause T e to fall below T g in the E-region, but Banks (1966a) pointed 
out that there the high gas density limits the electronic thermal conductivity to 
a lower value than Spitzer's expression gives and that allowing for this eliminates 
the negative values of T e - T . Above 500 km the temperature rises steadily to 
a value which, in the equatorial plane, is calculated to be 3000-5 000°K. Tempera- 
tures of this order have been observed by Serbu and Maier (1966) at altitudes of 
several thousand km. 

Banks (1967b) has studied the effect of thermal conduction amongst the ions. 
Below 500-600 km, there is little effect, T being slightly enhanced by heat con- 
ducted down from above. At greater heights T i remains significantly lower 
than T . In one example where T = 2600°K, T = 1000°K, T is 2500°K without 
and 225 0°K with thermal conduction. 

(b) Behaviour of Ionospheric Temperature by Day 

Taking all of the factors described in the last section into account, it becomes 
possible to calculate model profiles of T e and T x which can be compared with 
observation. Dalgarno et al (1968) made calculations of this sort for comparison 
with incoherent scatter results for Millstone Hill for 1963 July (Figure 17), The 
agreement is generally good except that the heat input must be reduced by about 
60% at 200 km where the calculated temperatures are too high. Nagy et al (1969) 
also made calculations for Millstone Hill and obtain temperatures which were 
rather too low near 200 km suggesting that the discrepancy may be quite sensi- 
tive to the details of the calculation such as the assumed solar spectrum and 
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photo-ionisation cross-sections. Hermanand Chandra (1969) made a similar 
computation for Arecibo, but obtained T e values which were slightly too low 
above 170 km. Lejeune and Petit (1969) calculated Q values for St. Santin pro- 
files of T and T from 120-500 km, and found the observed values to be less 
& 1 

than the calculated, which they attributed to the effect of non-local heating which 
they believed to be important down to 200 km. Despite these relatively small 
discrepancies, there is little doubt that the solar radiation provides the majority 
of the energy input to the day ionosphere. 

Another respect in which the calculations of Dalgarno et al (1968) failed to 
match the observations was that the computed temperatures show a smaller 
height gradient in the topside ionosphere than is observed. In these calculations, 
non-local heating and thermal conduction was neglected. In the calculations of 
Nagy et al, not only were these effects included, but so also was the energy 
balance of the magnetosphere. This allows both the energy conducted downwards 
from the magnetosphere and the photo -electrons from the conjugate hemisphere 
to be allowed for. When these contributions are included, a gradient in T e of 
2°K/km throughout the topside ionosphere is obtained. The downward flux of 
photo-electron energy from the conjugate ionosphere is found to be 2.8 x 10 9 
e'V/cm 2 sec, whilst the heat conducted downwards from the magnetosphere (which 
by 14 h LST, the time for which the calculation is made, is assumed to have 
reached a steady state with all the heat deposited being conducted downwards) 
is found to be 8 x 10 9 eV/cm 2 sec. Thus for this rather high magnetic latitude 
the greater part of the escaping photo-electron energy is deposited m the magne- 
tosphere. The proportion will be smaller at lower latitudes because of the 
smaller electron content of the geomagnetic field tubes. 
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The importance of thermal conduction and non-local heating are very clearly 
demonstrated by comparing the Millstone Hill profiles in Figure 18 with those 
for Jicamarca. Since the dip latitude at Jicamarca is only 1°, the magnetic field 
lines are nearly horizontal and these effects /are suppressed there. T e then rises 
to a maximum near 230 km, falling back to T g at greater heights precisely as 
predicted originally by Hanson, in a theory which neglected these effects. 

Geisler and Bowhill (1965a) showed by scaling various parameters in the 
calculation of T g that at high values of the collisional cooling rate in the 200- 
300 km range, the monotonically increasing temperatures shown in Figure 18(a) 
might be replaced by a variation similar to those of Figure 18(b) even at middle 

i 

latitudes. Such a behaviour could indeed be„seen in the profiles/ obtained on some 
early rocket flights [Spencer et al (1962)]. Results both from Arecibo [Mahajan 
(1967a)] and St. Santin [Carru et al (1967b)] show that on days when n e at the 
F peak exceeds (5-6) x 10 s cc" 1 a minimum in T e develops near 250-300 km as 
a result of the increased collisional cooling. On days when n e is lower a mono- 
tonic increase in T is observed. 

e 

The general features of the diurnal variation of temperature are shown in 
Figure 19. Throughout the day the variation is slight, with a fall during the night 
hours towards T g . A striking feature, shown by all reported observations in the 
topside ionosphere [Bordeau et al (1964), Evans (1965a), Brace et al (1965), 
Sagalyn et al (1965), Willmore (1965), Carlson (1966), Farley et al (1967), Brace 
et al (1967)] , at least during some epochs, is that the temperature maximum 
does not occur near or after noon, but at around 09 h local time. The magnitude 
of this early morning overshoot varies from 100-2000°K. The effect is greatest 
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between 300<and 500 km altitude and also near the magnetic equator where the 
overshoot is between 700 and 2000°K, falling at magnetic latitudes of 30°-50° to 
200-600°K. At high latitudes the effect may be confined to the summer months 
[Evans (1965a)] . This-dawn maximum appears to be the result of a more rapid 
build up of neutral gas density (and so of heat input) than of ionisation (or colli - 
sional cooling) during the early morning, so that the cooling rate lags relative to 
the heating rate. The inhibition of upward diffusion of ionisation near the mag- 
netic equator would then account for the enhanced effect seen there. That such 
an effect should occur was first predicted by Dalgarno and McElroy (1965); 
da Rosa (1966) solved the time-dependent heat balance equations for the dawn 
period and obtained results very similar to those which are observed. 

The principal geographical feature of the day electron temperature is an 
increase from equator to latitude 50° of about 600°K, after which the temperature 
slightly decreases (Figure 20) [Bowen et al (1964c), Brace et al (1967)] . The 
distribution shows marked geomagnetic control, being symmetrical about the 
magnetic not the geographic equator. Dalgarno suggested that this did not indi- 
cate an additional high latitude heat source but arose only from the decrease in 
n e and so of the collisional cooling rate. To this effect, at high altitudes, must 
be added that of the suppression of non-local heating near the magnetic equator. 
From an analysis of Ariel I data, Willmore (1964) showed this explanation to be 
correct. Superimposed on the latitude effect are smaller scale variations of 
±2Q0°K [Willmore (1965)]. Since these are invariably negatively correlated with 
n e variations that also occur they appear to arise merely from local anomalies 
in n e which produce variations in the collisional cooling rate. Thus all the geo- 
graphical features in the T e distribution can be explained by variations in elec- 
tron number density. 
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The behaviour of T i is also shown in Figure 18. It will be seen that T x rises 
from T g to T g as described in the previous section. In particular, values of T. 
greater than T e have not been observed, thus confirming the absence in general 
of a strong direct heat source for the ions. In the transition region from H + to 0 + 
the enhanced collisional cooling may produce negative temperature gradients 
[Banks (1967b)] which may possibly have been observed [ Wrenn and Smith (1969), 
Colin et al (1969)3 • - 

(c) Behaviour of Ionospheric Temperature By Night 

On the basis of the account presented in Section 6(a), T e would be expected 
after sunset to fall rapidly to T g . Whilst this is observed to happen at altitudes 
near 200 km, it is clear from Figure 20 that it does not do so at greater heights, 
but falls slowly to a minimum at 02-03 h . There is a good deal of evidence that 
T e > T g at many places through the night. It is found (Figure 18) [Bowen et al 
(1964c), Brace et al (1967)] that there is still a latitude increase at night although 
the temperature m inimu m over the equator is much broader than it is by day, 
the rise occurring at a geomagnetic latitude of about 30°. Considerable altitude 
gradients are observed from 400-800 km, especially in the winter when values 
of 0.9°K/km or more are recorded, indicating a significant downward heat trans- 
port [Prasad (1968), Mahajan and Brace (1969)] . Finally there are substantial 
seasonal and storm-time variations of T e [Evans (1965b), (1965a), (1967b)] , 
which do not result from variations in T alone. Thus, T > T at least for 

g eg 

some latitudes, altitudes and seasons and we are forced to admit the existence 
of a night time source of energy. Gei'sler and Bowhill suggested that this 
might be the energy which had been stored in the magnetosphere during the day 
when it is raised to a high temperature. Nagy et al (1968) using Geisler and 
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Bowhill's estimates of magnetospneric temperature solved the time-dependent 
heat balance equation and showed that hear fluxes exceeding 10 8 eV/cm 2 sec 
could be maintained for several hours after sunset. Nisbet (1968) predicted rather 
higher temperatures in the magnetosphere in agreement with observations [ Serbu 
and Maier (1966), Norman (1969)] so that Nagy's values may be an under-estimate. 
Moor croft (1969) determined the night time temperature gradients over Arecibo 
and deduced downward heat fluxes 85 min after sunset of 4.3 x 10 8 eV/cm 2 sec, 
declining 5 hours after sunset to 3.4 x 10 7 eV/cm 2 sec. It thus appears that 
magnetospheric heat storage is adequate at least to explain the night enhancement 
of T e , the latitude effect arising from the greater storage on longer field lines. 
Mahajan and Brace (1969) deduced the temperature gradients at night from near- 
simultaneous passes of Explorer XXII at 1000 km with Alouette II at 2500 km. 

The temperature gradients at 1000 km were deduced to be 0.6°K/km and 1.2°K/km 
for the 30° and 50° field lines, values which were considered to be consistent with 
magnetospheric heat storage. In addition if the conjugate ionosphere is illumi- 
nated, photo-electrons generated there may contribute to heating in the night- 
time, an effect which will be further discussed below. 

It has also been suggested [Willmore (1964), Spencer et al (1965), Brace 
et al (1965)] that precipitated or partly trapped electrons with kinetic energies 
of the order of 100-2000 eV might be responsible for the maintenance of night- 
time temperatures. Whilst magnetospheric heat storage is certainly the pre- 
dominant process, the possibility still remains that soft electron heating may 
produce additional heating at some times or places. Thus, Knudsen (1968) ob- 
served various regions of intense electron precipitation, the strongest being one 
near the S. Atlantic magnetic anomaly where the particle flux exceeded 10 11 /em 2 sec 


58 



and the energy flux 1 erg/ cm 2 sec. In this case it was concluded that the soft 
electrons were secondaries resulting from interaction of the radiation belt 
particles with the atmosphere. Other high flux regions coincided with areas 
where enhanced electron temperatures were observed in Ariel I. Probably these 
fluxes were unusually large since they much exceed those reported by Maier 
(1969). 

Dalgarno (1964) remarked that the excitation of the N t bands in the night 
airglow places an upper limit on the precipitated electron flux of 10" 2 ergs/cm 2 sec. 
Nathan and Seaton (1966) showed that at 400 km a trapped flux 100 times as great 
can exist, so that the heating may be due to an electron flux of 1 erg/cm 2 sec 
without undue airglow intensity. These may have a heating effect comparable 
with or greater than that of conduction from the magnetosphere and so would 
produce a significant effect on the night ionosphere. 

After the minimum in the early morning hours, the electron temperature 
starts to rise again somewhat before ground sunrise. During the summer months 
this rise occurs at sunrise at 100 km but in the winter an initial increase may be 
observed several hours earlier, followed by an increase in the heating rate at 
local dawn [Carlson (1966), Carru et al (1967b), Evans (1967b) - ] . The effect has 
been observed even at altitudes of 2500 km [Shepherd and Wrenn (1969)] . The 
time of the initial rise in T e correlates well with the time of ionospheric sunrise 
at the magnetically conjugate point which is of course in the summer hemisphere 
where sunrise is earlier. The effect is therefore ascribed to the production of 
photo -electrons in the summer hemisphere which, travelling through the mag- 
netosphere, produce heating in the winter hemisphere. This process was first 
suggested by Cole (1965) in connection with the predawn enhancement of the 
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X6300 airglow. In their observations of photo-electrons from the plasma line at 
Areeibo referred to previously, Yngvesson and Perkins (1968) indeed found down- 
ward fluxes of photo-electrons during this period. The pre-dawn heating rate 
will be a measure of the photo-electron flux; it ranges from 4-6°K/min at Areeibo 
and St. Santin down to l°K/min at Millstone Hill. The latter station is at a high 
magnetic latitude (L = 3.2) suggesting that a large proportion of the photo- 
electrons is absorbed in the long magnetospheric path for this station. This is 
in agreement with the calculations of Nagy et al (1969) for this location which 
showed 74% of the photo-electron energy to be absorbed in the magnetosphere. 

The coupling of the two hemispheres which results from the existence of 
the photo-electrons has a significant influence on the seasonal variation of the 
night temperatures which may become large at high latitudes [see Figure 11(b)]. 
The winter temperatures at Millstone Hill exceed the summer values by about 
800°K. Evans (1967b) has pointed out that the conjugate location to Millstone Hill 
is at 72°S, 81°W where the atmosphere remains sunlit throughout the day for 
several weeks during the northern hemisphere winter. This will result in high 
night-time temperatures in the magnetosphere and in sustained photo-electron 
heating which are presumably responsible for the high winter temperatures at 
Millstone Hill. Prasad (1968) has noted a similar seasonal variation in the night 
temperature gradient above 400 km at Areeibo. The temperature gradient is a 
maximum in winter, presumably because due to the late sunset in the summer 
hemisphere the magnetosphere has less time to cool. To this effect must be 
added the greater photo-electron heating in the winter hemisphere. 

It may be noted that these two effects, the maintenance of high electron 
temperatures at night and the predawn increase of T e , depend on a fortuitous 
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circumstance, since the first requires that the photo -electrons deposit a signifi- 
cant fraction of their energy in the magnetosphere whilst the second requires 
that the fraction should not be too large so that considerable energy still impinges 
on the conjugate hemisphere. It is very remar liable that magnetospheric densi- 
ties should just be such as to satisfy these opposing conditions over a wide range 
of latitudes and other conditions. 

8. IONOSPHERIC TEMPERATURES BELOW 140 km 

It will immediately be apparent from Figures 13(b) and 16 that below 140 km, 
Q is falling whilst the loss rate is increasing so that T e will be expected to ap- 
proach T We have also noted that the reduction of electronic thermal conduc- 
tion by the neutral gas prevents the ingress of heat from above. Before dis- 
cussing the observations, which unfortunately do not lead to an unambiguous 
picture, we will discuss some possible alternative mechanisms for maintaining 
enhanced electron temperatures. 

Walker (1968) has suggested that in the E-region the electrons will be in 
thermal equilibrium with the vibrational levels of N 2 . The vibrational tempera- 
ture of N 2 may be enhanced through eollisional deactivation of O^D) by N 2 . A 
large amount of solar energy is absorbed in the E-region in the dissociation of 
0 2 through the Schumann-Runge bands and the production of 0( 1 D) . Thus the 
total energy gained by N 2 is considerable but the fraction of it which goes to 
vibrational excitation is not known, Walker estimates that 1/3 - 1/7 of the total 
would suffice to maintain (T e - T ) at several hundred °K, though as fine structure 
cooling was neglected, this is an underestimate. Measurements of N 2 vibrational 
temperature can be made in aurorae when it is not invariably enhanced above T g , 
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so an additional quenching mechanism would be required on some occasions. T e 
is not invariably enhanced in the ionosphere, so a variable quenching mechanism 
would be required here also. On the other hand Dalgarno et al (1968) suggest 
that other collisional processes will always maintain the population of the vibra- 
tional levels at that corresponding to the gas temperature, so preventing, such an 
enhancement of T . 

e 

Dalgarno and Henry (1965) showed that in the D-region, T e would similarly 
be determined by the rotational temperature of N 2 which is presumably equal 
to T g . Photodetaehment from O' which produces electrons with a kinetic energy 
of 2.9 eV, a process originally proposed by Sears (1963), is not effective in 
raising T e . Moreover, complex ions based on NO 3 “ and not 0 2 ~ are the dominant 
negative ions of the D region. 

The effect of electric fields in the E-region has been discussed by Rees and 
Walker (1968). Their results will be discussed in more detail in a later section, 
but we may note here that electric fields of the order associated with the Sq 
current systems in the ionosphere, i.e. about 1 mV/m, will not produce a sig- 
nificant effect except in the vicinity of the equatorial electrojet [Satya Prakash 
et al (1968)] . At middle and high latitudes the principal effect of the electric 
field (assumed to be perpendicular to the magnetic field) is to raise the ion 
temperature. Near the equator the effect is strongly enhanced by the production 
of a vertical polarisation field whose principal effect is to heat the electrons. 

For an applied field of 1 mV /m and an electro jet current density of 5 amp/km 2 , 
the energy input is sufficient to raise T g by 2-300°K. Walker (1966) has also dis- 
cussed the D-region heating by electric fields of about 2.5 mV/m arising from 
neutral winds and concludes that the temperature rise here also will be insignificant. 
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The effects of wind motions in the E region has been studied by Gleeson 
and Axford (196V) who find that a windshear pattern which produces a strong E s 
layer with a peak density .of 8.5 times the ambient value would lead to a depres- 
sion of about 70°K in T e at the centre of the layer, a slight enhancement occurring 
above and below. Hooke (1969b) points out that the wind shears are themselves 
the result of atmospheric wave phenomena which will be accompanied by varia- 
tions in T g estimated to be ±50°K. When this is taken into account, it is found 
that T e in the E g layer may either be enhanced by 60°K or depressed by as much 
as 60°K relative to its value in the undisturbed ionosphere above and below the 
layer. 

In general, then, there is no clear indication on theoretical grounds of any 
mechanism which will maintain D and E region electron temperatures signifi- 

t 

cantly (i.e. 400°K) above the gas temperature, although temperature variations 
of a detectable magnitude may occur in the vicinity of the equatorial electrojet, 

E s layers or (as will be seen in Section 10b) in aurorae. 

No direct measurements of T have been made much below 100 km where 

e 

the occurrence of collisions vitiates the usual theory of the probe. JBetween 
100 and 140 km, the experimental evidence is, at least to some degree, conflicting. 
Figure 21 shows temperatures obtained during the day from three incoherent 
scatter stations [Carru et al (1967a), Evans (1967a), Ward and Perkins (1968)] . 

A complex altitude and diurnal variation is observed and attributed to the varia- 
tions of T g associated with internal gravity waves. For St. Santin the range at 
each altitude during daylight hours is indicated by a vertical line which thus indi- 

cates real variations as well as random errors which are smaller. It will be 

\ 

seen that T g obtained in this way and T g are in good general agreement up to 
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130 km. It is possible that T e is systematically higher than T g but if this is so, 
the difference does not exceed 50-100°K. 

Many rocket measurements of T q in the E region by Langmuir probes have 
been reported. Generally speaking, in the region from 100 to 140 km, values 
between 500°K and 1000°K are obtained. There is often a pronounced scatter of 
the measurements. Smith et al (1968) have summarised the results of 29 rocket 
flights at 145 km. They conclude that T e shows a general trend towards T g at 
the minimum of the solar cycle and that it shows considerable variability at any 
site with a time scale greater than several hours, but less than several days. 
There is thus a clear conflict with the incoherent scatter results in contrast to 
the situation at greater altitudes. The cylindrical probe comparison at Arecibo 
[Brace et al (1969)] is the only one where results at these altitudes were obtained, 
the comparison having been made down to 130 km without any indication of a dis- 
crepancy. Both methods are subject to the effect of collisions in the neutral gas 
at low altitudes. This necessitates a correction to the results of incoherent 
scatter below 120 km, and it is possible that this is not made with sufficient 
accuracy. In the case of the Langmuir probes, no theory which would allow a 
correction to be made is available and in fact the altitude at which a significant 
error may occur in T e cannot be identified with certainty. However, collisions 
become important for incoherent scatter when the electronic mean free path is 
comparable with the radio wavelength, and for the probe when it is comparable 
with the probe dimensions which are usually much smaller. Thus collisions 
should not usually be important for the probe until 100 km or less is reached 
and some versions of the probe should be much less affected than others. Thus 
the cylindrical probe, on account of its small diameter, should function to much 
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lower altitudes than the nose-tip probe. Yet in the observations discussed by- 
Smith et al (1968) which were obtained by both these methods, no such difference 
is apparent. Another possible source of error may be neglect of the geomagnetic 
field if for some reason its effect were more pronounced in the E-region. Once 
again, because of the relative sizes, the nose-tip probe should be more affected 
than the cylindrical probe. If the difference between probe and incoherent scatter 
arises from a- defect m one technique or the other of this land, it should always 
be present. In Figure 21 are shown probe results for four rocket flights, two 
each by the cylindrical probe and the ac probe which were shown in Section 6(b) 
to be in good agreement with incoherent scatter at greater heights. To take ac- 
count of the possibility of an error due to aerodynamic heating as suggested by 
Smith [Section 3(g)] , results from the downleg only of the trajectory have been 
used, although this is not obviously relevant for the probe systems employed 
here, and the ascent values were actually lower for two of the rocket flights. 

The ac probe results [ Aubry et al (1966)] were obtained during the afternoon at 
Hammaguir at almost the same local time two days apart. The rocket installa- 
tions were similar. On the first day, T e was found to be close to the CIRA 1965 
values of T g . They also agree well with the incoherent scatter values, tending 
to be on the low side of these. On the second flight two days later, T e values 
more than twice as great were obtained. Similar results were obtained from the 
second pair of flights [Brace et al (1969a)] . These were made from Fort 
Churchill. The first, during the day, showed temperatures systematically about 
100°K above T , .whilst the second which occurred during the subsequent night 
showed results nearly twice as great. The high temperatures on this occasion 
may easily be explained in terms of excitation processes in the auroral zone 
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[see Section 10(b)] ; the day results show that they are certainly not the result of 
a systematic error in the cylindrical probe technique. 

Since there are certainly some occasions when relatively low temperatures 
are observed by Langmuir probes, a defect in the theory of operation which 
would lead to a systematic error can be ruled out. Moreover, in the cases cited, 
there appears to be no reason to believe that the measurements leading to high T e 
values were less satisfactory in quality than the others. In rocket measurements 
of T e m general, enhanced values are more common than values close to T , 
though it is often noticeable that there is a large scatter in T e , with the lowest 
values close to T g , rather than being symmetrically distributed about T g as 
might occur simply from random experimental errors. Since high values of T e 
are apparently not observed at all by incoherent scatter, the discrepancy between 
the two methods can hardly be an accidental effect of sampling alone. Thus it 
seems that the difference is real and yet does not arise simply from an inade- 
quacy in the experimental technique in either case. A possible explanation is 
that T g enhancements which are localised occur, so that T e exhibits a very spikey 
variation. The high spatial resolution of the Langmuir probe would enable these 
to be partly or completely resolved. This would account for the large scatter 
often observed and for the fact that despite this scatter T e values lower than T g 
are not observed. The radar, however, has a relatively poor spatial resolution 
ranging from 3 km at St. Santin to 30 km at Millstone Hill and so may measure 
the mean value of these variations, which need be only moderately greater than 
T g . Alternatively it may be that in the E-region, despite theoretical predictions 
to the contrary, the electron energy distribution is non-Maxwellian. If tempera- 
ture is an unsuitable parameter to describe the energy distribution, it would not 
be surprising if different techniques gave discordant estimates of it. 
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If it is supposed that small scale enhancements of T e above T g occur (though 
it should be noted that high values of T e might, in principle at least, reflect 
similar small scale enhancements of T g , presumably such strong turbulence does 
not exist in this region) then a source of energy to maintain the difference must 
be found. We have seen that enhanced values were observed at night at Fort 
Churchill. Moreover, Smith et al (1965, 1968) made observations in two eclipses 
(1963 July 20 and 1966 Nov. 12) and found no evidence for a fall during the eclipse. 
Thus it appears that a direct effect of solar radiation is not responsible unless a 
molecule with a long life time excited state is involved. Walker (1968) states 
that the excited vibrational states of N 2 will persist throughout the night so that 
high temperatures might be maintained both at night and through an eclipse. 

A much smaller number of measurements of T has also been made between 

1 

100 and 140 km by both techniques with more consistent results. Knudsen and 
Sharp (1965) using an ion trap on a rocket fired near local noon, observed a wave- 
like variation of T i between 100 and 160 1cm, the amplitude of the variation being . 
100°K and the wavelength 12 km. The mean value of T x was only 100°K above the 
expected T . They suggested that the variation of reflected the variations 
of T g resulting from density variations in intense gravity waves, although waves 
such short wave-length would be strongly damped by thermal conduction, and the 
wave amplitudes also were an order of magnitude greater than Hines (1965) had 
estimated for such waves in the E-region. Using a similar method, Harris et al 
(1969) observed wavelike variations of T from a satellite near 200 km altitude. 
These were not invariably present and were generally observed at high latitudes. 
The wavelength was around 200 km and the wave amplitude 13% of the mean tem- 
perature. Once again the wavelengths are shorter than expected for gravity waves. 
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Wand and Perkins (1968) using incoherent scatter at Arecibo observed the 
variation of T (presumed equal to T e ) from 97-128 km. The results showed an 
oscillation of peak amplitude 50°K about the CERA values of T . Hooke (1969a) 
suggested that these might be associated either with internal gravity waves or 
with tidal motions of the neutral atmosphere, and that only the latter was con- 
sistent with the temporal stability of the observed patterns. 

Vasseur and Waldteufel (1969) observed both plasma drift motions and T r at 
St. Santin. They showed that these were consistent with the occurrence of in- 
ternal gravity waves with a period of 2 h. The peak amplitude of the variation 
of T 1 was 50°K. The vertical wavelength was found to be 520 km. Correlated 
variations in N e at the F2 peak were observed, a fact which was utilised to study 
the horizontal variations using ionosonde data. Thus it was deter min ed that the 
horizontal wavelength was 4000 km and that the waves propagated southward 
across Europe from the auroral zone where it was suggested that they were 
generated. 

Thus there appears to be quite good evidence for the observation of tempera- 
ture variations of ±50°K or a little more, associated with atmospheric movements 
in internal gravity waves and possibly also of both shorter and longer scales than 
these waves. It also appears that the gravity wave pattern has a major influence 
on the apparent diurnal variation of temperature in the E region [Carru et al 
(1967a)] . 

9. SOLAR CYCLE VARIATIONS IN E-REGION IONOSPHERIC TEMPERATURE 

The calculations of Geisler and Bowhill (1965a) suggested that the solar cycle 
variations in T e would not be very large. From the results shown in Figure 11 
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it is evident that, except under circumstances where T e is not too much greater 
than T n so that it shows a similar variation, the day-to-day variability from 
other causes such as ionospheric storms is such as to mask the solar cycle 
changes unless careful averaging of substantial numbers of observations is car- 
ried out. This will best be done with results from incoherent scatter stations 
when these have operated for a sufficient period. 

Attempts have been made to correlate changes m T g with variations m solar 
EUV flux over shorter periods using the 10.7 cm radio emission as a monitor of 
the latter. Willmore (1965) found a variation of (4.6 ± 2.3) °K/flux unit* in the 
day at 500 km. Here, the 10.7 cm flux on the day of the T e observation was used. 
Mahajan (1967b) correlated Arecibo values of T e for 1965-7 with 27 day averages 
of 10.7 cm emission, obtaining a night-time change of 3.4°K/flux unit and, at 
noon, below 400 km, 4.5°K/flux unit. Since it is found for variations in exospheric 
temperature that the coefficient relating changes during one solar rotation 
(1.9°K/flux unit) is less than that for long-term' variations (3.46°K/flux unit), a 
difference which presumably reflects the real behaviour of the EUV and radio 
emissions, the two results are not really in as good agreement as at first ap- 
pears. Nonetheless, the variations in T e appear to be similar to those in exo- - 
spheric temperature. 

10. IONOSPHERIC TEMPERATURES UNDER DISTURBED CONDITIONS 

(a) Ionospheric Storms at Moderate and Low Latitudes 

The behaviour of T e in several small ionospheric storms was studied using 
Ariel I observations by Willmore (1965), who found rather localised reductions 

*The flux unit is 10 -22 'W/cm 2 sec Hz. 
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to occur at many places on the globe, the decreases amounting to 100°K or more. 
The behaviour was much the same by day or night and the decreases- in T e were 
well correlated with increases in electron density. Only on one occasion, at 
4 locations near to the S. Atlantic magnetic anomaly was an increase in T e - re- 
corded. Knudsen (1968) made observations with a retarding potential analyser 
of T , n e and 10 eV electron flux. -He found intense particle fluxes (exceeding 
10 11 el/cm 2 sec) near the magnetic anomaly which he attributed to secondaries 
produced by particles of the radiation belt. Two of the high flux regions coincided 
with the region of enhancement observed in Ariel -I. Reddy et al (1967) obtained 
very similar results from Explorer XXII though the depressions in temperature 
which they observed were larger in magnitude, reaching 300°K. It seems that 
this behaviour is typical of rather weak storms, and possibly of the initial phase 
of more intense storms [Bauer and Krishnamurthy (1968)] . 

It was suggested by Willmoredhat the change in T e resulted largely or com- 
pletely from changes in the collisional cooling rate consequent on the enhancement 
of n e . Thus, the principal change in the storm is that of n e , T e changing there- 
after in a maimer dictated by the energy balance of the ionosphere. In particular, 
no additional energy source than that provided by solar radiation need be present. 
Reddy et al advanced an alternative explanation, that the photoelectron flux is 
reduced in regions of enhanced electron density because of increased scattering 
which raises the escape level. However, if T e and n e both increase, as in the 
regions near the magnetic anomaly, it is at least possible that an additional 
source of energy must be postulated in these regions. 

Evans (1965) made observations through two storms of moderate severity 
at Millstone Hill in 1963 September. T x was observed to be little affected except 
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below 400 Ion where small increases of 100-200°K occurred. T however was 

e 

enhanced substantially, by 700-800°K. At night this implied an increase in the 
energy input by a factor of three. At Millstone Hill, where most of the photo- 
electron energy is normally deposited in the magnetosphere, it is doubtful 
whether the magneto spheric heat output can be increased by so large a factor. 

Henderson and Wrenn (1970) used Explorer XXXI observations to study 
several storms early in 1967. They found that both enhancements and depres- 
sions of considerable magnitude, i.e. 500-1000°K, could occur at all heights from 
500-1600 km at least. In a given storm, different effects may be observed at 
different locations. Thus in the storm of 1967 June 6 T e increased by 400°K 
over S. America together with a large increase of electron density. A few orbits 
later over Africa it was found that a similar increase of T had occurred there 

e 

but with only minor decrease in density. Such cases, where T e increases and n e 
either increased or only decreases slightly, are suggestive of an enhancement 
in the energy input and so of an additional energy source besides solar radiation. 
However, Henderson and Wrenn also observed very substantial changes in ion 
composition during storms, which will result in changes in the collisional cooling 
rate. Atmospheric temperature changes must also be allowed for. It follows 
that very careful calculations are needed before it is possible to decide what - 
changes have actually occurred in the energy input. No clear pattern in the 
temperature changes in ionospheric storms emerges from the results available 
so far. Because these changes, which will result in redistribution of ionisation 
along the geomagnetic field lines and may also be accompanied by changes in the 
loss rate near the F2 peak, can considerably modify the morphology of the 
ionosphere, it is desirable to obtain a clear picture of them as a contribution to 
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building a theory of the ionospheric storm. Moreover-, if an additional source of 
ionospheric energy exists during the storm, it would be most important to iden- 
tify it. The question of the existence of such a source seems one of the most 
significant unanswered questions concerning ionospheric temperatures in the 
T -region. 

(b) Ionospheric Temperatures in the Auroral Zone 

In the neighbourhood of an aurora and in the auroral zone generally addi- 
tional heating processes are likely to modify substantially the picture previously 
given, since large amounts of energy are certainly available in the production of 
aurorae and the associated high latitude current system and some part of this is 
likely to be dissipated as ionospheric heating. The energy loss mechanisms 
under these circumstances will not differ from those that have already been 
discussed so that it is only necessary to describe the heating processes. 

The first of these to be proposed was Joule heating from electric fields 
produced in the magnetosphere [Cole (1962), Megill et al (1963), Walker (1966), 
Cummings and Dessler (1967), Rees and Walker (1968)3 • Considerations of the 
magnetic field changes in storms show that the electromagnetic energy density 
available over the auroral zone is not far short of that due to solar EUV radia- 
tion in the sunlit hemisphere, so this is potentially a significant source. Com- 
ponents of the electric field E perpendicular to the magnetic field have been 
measured by barium cloud releases over Kiruna [Foppl et al (1968)] , values in 
the range 2-20 mV/m being obtained. A rather larger value of 50 mV/m is re- 
quired to account for the current flowing in the auroral electro] et [Bostrom 
(1964)] . At high latitudes where the dip angle is large, the polarisation field 
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which, occurs in the equatorial electro] et is absent and the heating process is 
rather different from that which occurs near the magnetic equator- In perpen- 
dicular electric and magnetic fields, similar drift velocities are imparted to 
ions and electrons, so that unlike all the other mechanisms discussed the energy 
gain of the ions in this process is much larger than that of the electrons [Rees 
and Walker (1968)] . Also, the heating effect is evidently greatest when the ion 
collision and gyrofrequencies are comparable, i.e. in the E-region. In conse- 
quence of these two facts, T x rises above T e , this rise commencing in the 
E-region and being maintained into the F-region. T g is in fact largely enhanced 
not as 'a result of the direct heating by the electric field, but as a result of ion- 
electron collisions just as the ion temperature is enhanced in the normal F- 
region. For electric fields up to 50 mV/m; the effect on T e is small, but T 
rises to 2000°K over the range 160-300 km. 

An obvious source of heating under auroral conditions is the precipitated 
electron flux responsible for the optical excitation. Cole (1967) showed that such 
a precipitated flux would generate a space charge resulting in the flow of field 
aligned currents, and that the accompanying Joule heating would exceed that of 
the primary electrons in some circumstances. Walker and Rees (1968a) and 
Stolarski' (1968) pointed out that the incident beam of primary particles in the 
atmosphere is accompanied by a shower of low energy secondaries produced by 

i 

ionisation of the neutral gas. Because of the increase of the Coulomb and other 
relevant cross-sections at low energies, the heating produced by the secondary 
electrons is very effective and, in fact they and not the field-aligned currents 
represent the dominant energy source. They clearly play the same role as the 
photoelectrons in the heating by solar EUV photons, their initial energy being 
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degraded by collision in a similar way. They have a heating efficiency ranging 
from a fraction of an eV at low altitudes to several eV at 400 km. Rather than 
starting directly from an assumed flux and energy distribution of precipitated 
electrons. Walker and Rees considered five auroral arcs ranging m brightness 
from 7.2 to 34.7 kR, for each of which a luminosity profile of the 3914A band 
had been measured. Since the optical luminosity enables the excitation rate to 
be determined, and since the relation between the excitation and ionisation cross- 
sections is known, the volume ionisation rate can be determined directly from 
the luminosity profile. This has only to be multiplied by the value of the heating 
efficiency for the appropriate altitude to obtain the heat input. Then T G and T 1 
can be calculated in 'the usual way (Figure 22). Maximum temperatures of 2900- 
3500°K were found to occur at 300-400 km. In these calculations, non-local 
heating was not considered, the temperature above 400 km being maintained 
largely by thermal conduction. The inclusion of non-local heating would pre- 
sumably remove the temperature maximum in the F-region, resulting in the 
maintenance of high temperatures throughout the topside ionosphere. More 
seriously, cooling by excitation of the fine structure levels of atomic oxygen was 
also neglected, so that the temperature rise would be seriously overestimated 
at low altitudes, though the error in the maximum value might not be so great. 

An important feature of the process is that as the precipitated electron flux in- 
creases both the energy input and the collisional cooling rate, which depends on 
n e and hence the ionisation rate, increase; hence the temperatures are not very 
sensitive to the precipated flux and will be much the same in intense as in weak 
aurorae. 
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Thus we see that the preeipated electrons responsible for the production of 
aurorae will raise T e rather than T ±> producing an enhancement to 2900~3500°K 
at 300-400 km. Higher values than this are unlikely to be produced simply by 
the particles even in very intense fluxes. If in addition an electric field E is 
present Tj will be raised with relatively little effect on T e . Since the energy 
input varies as E 2 , the enhancement will increase rapidly with E, reaching values 
of 2000°K at a fairly large value of 50 mV/m. The two effects 'are nearly inde- 
pendent, and measurements of T e and T x would enable them to be distinguished. 

Some measurements axe available which seem to confirm these predictions 
at least in part. Ulwick et al (1964) measured electron temperatures from a 
rocket through an auroral form over Fort Churchill, and found maximum values 
of 3500°K. Baker et al (1966) made measurements in a discrete auroral form, 
finding temperatures of 3500°K in the aurora and lower values of 2000°K outside. 
In a diffuse aurora, T e was no higher than 1800°K, the value rising with altitude 
up to the peak of the rocket trajectory at 193 km. Raitt (1969) made .observations 
with an ae probe on ESEO IA at an altitude of 800 km during the event of 1969 
Feb. 25-26, which was of moderate intensity. He found T e enhancements over 
the night time auroral zone of only 500°K. 

One major disturbance has been unusually well' observed, results having 
been reported from three satellites, Explorer XXH, Ariel HI and Explorer XXXI 
[Findlay et al (1969), Loftus (1968) and Henderson and Wrenn (1969)] . This was 
a large storm on 1967 May 25-26 when aurorae were visually reported from 
many locations in the US near an invariant latitude of 50°. The most remarkable 
observations were those of Findlay et al, who found that T e was enhanced from 
a pre-storm value of 2200°K to 6500°K near invariant latitude 50 D . In the night 
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hours, where reports of the position of the aurora were available, it was con- 
cluded that the temperature enhancements were directly above the aurora. The 
electron number density was also enhanced by a factor of 2. In the day, the 
enhancement of T e was less striking and n e was depressed. At invariant lati- 
tudes greater than 65° the enhancement of T. was much smaller and n e was 
again depressed. Considerable fine structure in T e was observed. Henderson 
and Wrenn detected no storm effect at all over Alaska, in agreement with these 
results. 

It is quite clear that the theoretical discussions described thus far would 
not account for the occurrence of temperatures as high as 6500°K. Whilst it is 
not certain that the event observed by Findlay et al is of this type, the most 
relevant study seems to be that of the mechanism of the stable auroral red are 
(SACARC) by Walker and Rees (1968). The principal feature of these arcs is the 
strong emission of the 6300/6364 A lines relative to other wavelengths thus re- 
quiring preferential exitation of the 0( 1 D) level. It is suggested that a high elec- 
tron temperature is responsible, excitation arising from electron impact in the 
high energy tail of the Maxwellian velocity distribution. This will give strong 
excitation of 0( i D) on account of its relatively low excitation energy. Walker 
and Rees discussed the following heating mechanisms: 

(a) an electric field E 

(b) heat conducted from high temperature regions in the magnetosphere 
and (c) precipitated electrons. 

They attempted to reproduce the characteristics of a particular red arc whose 
luminosity profile was recorded at Fritz peak, Colorado on 1961 June 22. It 
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was found (Figure 23) that an electric field of 180 mV/m would be required, 
being raised to 12000°K at 200 km, falling by 700 km to equal T e , which would 
reach 5000°K from 200 to 1000 km. However, a very large amount of energy 
would be imparted to the neutral atmosphere which would be raised to 3400°K in 
the arc, a value which Walker and Rees considered improbable. Moreover, such 
large electric fields have not yet been observed. Heat conducted from the mag- 
netosphere in a flux of (4-10) x 10 10 eV/em 2 sec or about an order of magnitude 
more than that due to the photoelectrons, would also suffice, raising T g at 
1000 km to 6000°K. The existence of magnetospheric energy source to maintain 
this supply of heat is of course hypothetical. If precipitated electrons are re- ‘ 
sponsible, they must be of such low energy that direct excitation of the line 
emission is small. Thus, electrons of about 15 eV in energy in fluxes of more 
than 10 9 /em 2 sec would be required. Walker and Rees note that the distinction 
between such low energy electrons and the heat conduction mechanism is not 
altogether clear. However, Maier (1969) has observed soft electron fluxes of 
the required magnitude over the auroral zone. 

Thus, if it is accepted that the electric field hypothesis is untenable, either 
thermal conduction of magnetospheric energy or a precipitation of very soft 
electrons from the magnetosphere are at least capable of reproducing the very 
high temperatures observed by Findlay et al. Measurements of T would show 
whether or not electric fields were important in this case. 

Rees and Walker (1968) also noted that even in the absence of visible aurora 

/ 

the zenith intensity of the 5577 A airglow is usually enhanced relative to mid- 
latitude values by a factor of 10, and suggested the existence of a steady excita- 
tion source' at all times. This may tend to raise T e and T x and certainly Brace 
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et al (1969a) and Nagy and Walker (1967) observed rather high night temperatures 
over Ft. Churchill. Rees and Walker found that a rather weak flux of precipitated 
electrons might give T e - T g rising from zero at 120 km to 600-800°K at 300 km, 
whilst an electric field of 50 mY/m would give - T g rising from 400°K at 
120 km to 850°K at 300 km, the two effects being nearly independent. T e enhance- 
ments of this magnitude are in good agreement with the rocket measurements, 
except that, as at lower latitudes, the latter usually show enhanced temperatures 
below 140 Ion which are apparently no easier to account for as the effect of 
precipitated electrons than as the effect of solar radiation. 
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FIGURE CAPTIONS 


Figure 1(a). Cylindrical probe characteristic obtained from ISIS-1 at latitude 26°, 
longitude 72°, altitude 600 km at 21 h 33 m GMT on 1969 October 4. 

Figure 1(b). Exponential portion of the characteristic after subtraction of the 

ion current, showing the determination of T e , i eo and space potential. 

• 

Figure 1(e). Actual and ideal characteristics in the accelerating rdgime, the 
ideal characteristic being calculated using T e and i eo determined from (6). 
[In practice, current ranges of different sensitivity would be used for the 
temperature and density determinations, leading to a better resolution of 
each quantity.] (Courtesy of Mr. L. H. Brace.) 

Figure 2. Explorer XXXI, showing the Langmuir probes. A spherical ion probe 
is mounted at the top on a short boom. On the left front face is the ac 
electron probe, on the right front face the screened plane probe; on the 
right side projects a cylindrical probe.* 

Figure 3. First and second derivatives of a spherical ion probe characteristic 
for an ion mixture, obtained on Explorer XXXI. [From Wrenn (1969)] 

Figure 4. Comparison of measured and theoretical second derivative charac- 
teristic for an 0 + plasma. The rms deviation of the measured points is 1.4% 
of the peak amplitude. [From Wrenn (1969)] 


Figure 5. The wake of a satellite (Ariel I) observed close to the surface. The 
relative currents to identical probes mounted on the base of the spacecraft 
and on a boom in relatively undisturbed plasma are shown as a function of 



the attack angle, i.e. the angle between the spacecraft axis and the velocity 
vector. Because the probe is in the base, it is on the forward or ram side 
of the spacecraft when the angle of attack is 180° and in the wake when this 
angle is 0°. The predominant ion was 0 + . [From Samir and Willmore (1965)] 

Figure 6. Ratio of electron densities measured with a probe in the surface of 
Explorer XXXI and those obtained by the topside sounder on Alouette II. The 
low densities measured at large Debye lengths are the result of the limitation 
of the probe current by a potential barrier caused by the proximity of the 
spacecraft. At about X D = 0.8 cm the effect disappears when a sheath forms 
about the spacecraft. [Courtesy of Dr. G. W. Wrenn.j 

Figure 7. Normalised incoherent scatter spectra calculated for an 0 + plasma 
(K = (M/m) 1/5 = 172) for propagation parallel to the magnetic field or no 
magnetic field. For X p = o.l (kX D = o.l) the spectral distribution is little 
affected by the presence of the ions. For X p = 3.0, the plasma line can be - 
seen at X = 3. [FromHagfors (1962)] 

Figure 8. Computed incoherent scatter spectra for ratios of the parameter T e /T i 
from 1.0 to 4.0. A plasma frequency of 5 MHz and' an ion temperature of 
1500°K were assumed. The results apply for 0 + ions and were made for a 
radar carrier frequency of 440 MHz and a pulse length of 1 msec. [ From 
Evans (1964)] 

Figure 9. Computed incoherent scatter spectrum for propagation almost per- 
pendicular to the magnetic field. The angles are (top) 87° and (bottom) 88°; 

I 

the dotted line shows the spectral distribution normal to the field. The 
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calculations are for a radar frequency of 40 MHz and typical ionospheric 
conditions. [From Farley, Dougherty and Barron (1961)'] 

Figure 10. T e and T 1 measured by the ac probes on Explorer XXXI compared 
with simultaneous and near-coincident observations by the St. Santin inco- 
herent scatter station on 1968 April 23. [From Henderson, Wrenn and 
Waldteufel (1969)] 

v 

Figure 11. Exospheric and electron temperatures determined by satellites and 
by incoherent scatter. The results are all for an altitude of 400 km and are 

for latitudes near 0° and 40° geomagnetic. The bracketed Explorer XVII 

/ 

points show the results of adjustment according to Eq. (23). The large vari- 
ation of the night temperatures at Millstone Hill is a seasonal effect." The 
vertical lines on the Arecibo observations indicate the range encountered 
and are not an error indication. [Data for these figures was obtained from: 
Explorer I — Jacchia and Slowey (1966); Jieamarca — Farley, McClure, 
Sterling and Green (1967); Ariel I — Willmore (1965); Explorer XVH — 
Brace, Spencer and Dalgarno (1965); Millstone Hill — Evans (1965a, 1967b); 
Arecibo — Mahajan (1967a), Prasad (1968).] 

Figure 12. (a) Variation of gas temperature T g with altitude in the ionosphere. 

(b) Variation of exosphere temperature with geographic position. 
[From Jacchia (1965)] 

Figure 13. Variation of (a) the electron heating efficiency e and (b) the.heat 
input Q with altitude at noon. (Non-local heating is neglected.) . [ From 
Dalgarno, McElroy,*Rees and Walker (1968)] 



Figure 14. Photoelectron energy distribution for various altitudes from 200- 

t 

300 km. [From Hoegy, Fourrier and Fontheim (1965)] 

Figure 15. Variation with altitude of the heat input Q showing both the local 

heating and also the additional non-local heating which results from including 
the contribution from photoelectrons originating below 300 km. [ From 
Geisler and Bowhill (1965)] 

Figure 16. Calculated electron cooling rates by various collision processes for 
1963 July at noon. The curves refer to excitation of (a) vibrational levels 
of N 2 , (b) rotational levels of N 2 , (c) fine-structure levels of O, (d) meta- 
stable levels of O, and (e) elastic collisions with positive ions. The cooling 
due to elastic collisions with neutrals and from rotational excitation of 0 2 
are negligible. [From Dalgarno, McElroy, Rees and Walker (1968)] 

Figure 17. Variation of electron temperature with altitude at Millstone Hill for 
1963 July at various local times. The solid curves are calculated including 
fine-structure cooling by atomic oxygen whilst it is neglected in the dash- 
dot curves. The dashed curves are values measured by Evans by incoherent 
scatter. [From Dalgarno, McElroy, Rees and Walker (1968)] 

Figure 18. Altitude profiles of T e and T x taken (a) at Jicamarea and (b) at Mill- 

v 

stone Hill. The effect on T e of the virtual suppression at Jicamarea of non- 
local heating is very striking. [From Farley, McClure, Sterling and Green 
* (1967) and Evans (1965)] 

Figure 19. Diurnal variation of T e at Arecibo in' winter. Also shown are times 
of sunrise and sunset at Arecibo and at the magnetically conjugate location. 
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Whilst at 225 km T e falls sharply at local sunset, and rises again at sunrise, 
at 400 km it declines steadily to a minimum in the early morning hours, 
starting to rise again at sunrise in the conjugate ionosphere. The early 
morning maximum- is most pronounced at 300 km on this occasion. [From 
Carlson (1966)] 

Figure 20. Latitude variation of T e as determined by probe measurements on 
Ariel I [ (a) and (b)] and Explorer XX [ (c) and (d)] . The Explorer XX results 
refer to an altitude of 1000 km. [From Willmore (1965) and Brace, Reddy 
and Mayr (1967)] 

Figure 21. Measurements of T e in the E-region by Langmuir probes and inco- 
herent scatter compared with gas temperature from CIRA 1965 (solid lines). 
6 and 7 refer to similar rocket flights, 6 near midnight and 7 near noon. 

8 and 9 were also similar, the two flights occurring at approximately the 
same local time, two days apart. 

Figure 22. Calculated temperature profiles for an auroral arc with \3914 A 

brightness of 34.7 kR. [T^ is the result neglecting thermal conduction, from 
Walker and Rees (1968a)] 

Figure 23. Calculated temperature profiles in a SAC ARC when the heating is 
supposed to result from (a) electric fields of various magnitudes normal to 
the magnetic field and (b) heat fluxes conducted from the magnetosphere. 

The temperature gradients are those at 1000 km. [From Walker and Rees 
(1968b)] 
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Figure 1(b). Exponential portion of the characteristic after subtraction of the ion current, 
showing the determination of T , i eo and space potential. (Courtesy of Mr. L. H. Brace) 
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Figure 1(c). Actual and ideal characteristics in the accelerating regime, the ideal characteristic 
being calculated using T e and i eo determined from (6). [In practice, current ranges of different 
sensitivity would be used for the temperature and density determinations, leading to a better 
resolution of each quantity.] (Courtesy of Mr. L. H. Brace) 






Figure 2. Explorer XXXI, showing the Langmuir probes. A spherical ion probe is mounted at the 
top on a short boom. On the left front face is the ac electron probe, on the right front face the 
screened plane probe; on the right side projects a cylindrical probe. 
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DECEMBER 16, 1965 10-31-57 Z 
ALT. 527 Km: LAT. 61°S: LONG. 68°W 




Figure 3. First and second derivatives of a spherical ion probe characteristic 
for an ion mixture, obtained on Explorer XXXI. [From Wrenn (1969)] 
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Figure 4. Comparison of measured and theoretical second derivative characteristic for 
an 0 + plasma. The rms deviation of the measured points, is 1 4 % of the peak amplitude. 

[From Wrenn (1969)] 
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i boom 



Figure 5. The wake of a satellite (Ariel I) observed close to the surface. The relative currents 
to identical probes mounted on the base of the spacecraft and on a boom is relatively undisturbed 
plasma are shown as a function of the attack-angle, i.e. the angle between the spacecraft axis 
and the velocity vector. Because the probe is in the base, it is on the forward or ram side of the 
spacecraft when the angle of attack is 180° and in the wake when this angie is 0°. The pre- 
dominant ion was 0 + . [From Samir and Will more (1965)] 
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Figure 6. ' Ratio of electron densities measured with a probe in the surface of Explorer XXXI and 
those obtained by the topside sounder on Alouette II. The low densities measured at large Debye 
lengths are the result of the limitation of the probe current by a potential barrier caused by the 
proximity of the spacecraft. At about \ D - 0.8 cm the effect disappears when a sheath forms 
about the spacecraft. ['Courtesy of Dr. G. W. Wrenn.] 





Figure 7. Normalized' incoherent scatter spectra calculated for an 0 + plasma (K = (M/m)^ 2 = 
172) for propagation parallel to the magnetic field or no magnetic field. For X - 0 1 (k/V D = 0.1) 
the spectral distribution is little affected hy the presence of the ions. For X** - 3.0, the plasma 
line can be seen at X = 3. [From Hagfors (1962)3 P 
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Figure 8. Computed incoherent scatter spectra for ratios of the parameter T e /T ( 
from 1 0 to 4.0. A plasma frequency of 5 MHz and an ion temperature of 1500°K 
were assumed. The results apply for 0 + ions and were made for a radar carrier 
frequency of 440 MHz and a pulse length of 1 msec. [From Evans (1964)] 









Figure 9. Computed Incoherent scatter spectrum for propagation almost perpendicular to the mag- 
netic field. The angles are (top) 87° and (bottom) 88°; the dotted line shows the spectral distri- 
bution normal to the field. The calculations are for a radar frequency of 40 MHz and typical 
ionospheric conditions. [From Farley, Dougherty and Barron (1961)] 




DEGREES K 

Figure 10. T e and T measured by the ac probes on Explorer XXX! compared with simultaneous 
and near-coincident observations by the St. Santin incoherent scatter station on 1968 April 23. 

[From Henderson, Wrenn and Waldteufel (1969)] 
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■ Figure 11. Exospheric and electron temperatures determined by satellites and by incoherent scatter. The results are all for an alti- 
tude of 400 km and are for latitudes near 0° and 40° geomagnetic. The bracketed Explorer XVII points show the results of ad|ustment 
according to Eq. (23). The large variation of the night temperatures at Millstone Hill is a seasonal effect. The vertical lines on the 
Arecibo observations indicate the range encountered and are not an error indication. [Data for these figures was obtained from. Ex- 
plorer l-Jacchia and Slowey (1966), Jicamarca - Farley, McClure, Sterling and Green (1967), Ariel I - Wtllmore (1965), Explorer XVII 
- Brace, Spencer and Dalgarno (1965), Millstone Hill - Evans (1965a, 1967b), Arecibo - Maha|an (1967a), Prasad (1968).] 
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Figure 13. Variation of (a) the electron heating efficiency e and 
(b) the heat input Q with altitude at noon. (Non-local heating is 
neglected.) _[From Dalgarno,, McElroy, Rees and Walker (1968)] 
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Figure 14. Photoelectron energy distribution for various altitudes from 200-300 km. 
[From Hoegy, Fourrier and Fontheim (1965)] 
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Figure 15. Variation with altitude of the heat input Q showing both the local heating and also the 
additional non-local heating which results from including the contribution from photoelectrons 
originating below 300 km. [From Geisler and Bowhill (1965)] 
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Figure 16. Calculated electron cooling rates by various collision processes for 1963 July ot noon. 
The curves refer to excitation of (a) vibrational level s«of N 2 , (b) rotational levejs of N 2 , (c) fine- 
structure levels of 0, (d) metastable levels of 0, and (e) elastic collisions with positive ions. 
The cooling due to elastic collisions with neutrals and from rotational excitation of 0^ are 
negligible. [From Dalgamo, McElroy, Rees and Walker (1968)] 
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MILLSTONE HILL, JULY 1993 



Figure 17. Variation of electron temperature with altitude at Mill- 
stone Hill for 1963 July at various local times. Tfie solid curves 
are calculated including fine-structure cooling by atomic oxygen 
whilst it is neglected in the dash-dot curves. The dashed curves 
are values measured by Evans by incoherent scatter. [From 
Dafgamo, McElroy, Rees and Walker (1968)3 
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Figure 18(a). Altitude profiles of T e and T ( taken at Jicamarca. The effect on T e of the virtual 
suppression at Jicamarca of non-local besting is very striking. [From Farley, McClure, Sterling 

and Green (1967) and Evans (1965)] 
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Figure 19. Diurnal variation of T e at Arecibo m winter. Also shown are times of sunrise and 
sunset at Arecibo and at the magnetically conjugate location. Whilst at 225 km T e falls sharply 
at local sunset, and rises again at sunrise, at 400 km it declines steadily to a minimum in the 
early morning hours, starting to rise again at sunrise in the conjugate ionosphere. The early 
morning maximum is most pronounced at 300 km on this occasion. [From Carlson (1966)] 
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Figure 20(a and b). Latitude variation of T e as determined by probe measure- 
ments on Ariel I. [From Willmore (1965) and Brace, Reddy and Mayr (1967)] 
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Figure 21. Measurements of T e in the E-region by Langmuir probes and incoherent scatter compared with gas 
temperature from CIRA 1965 (solid lines). 6 and 7 refer to similar rocket flights, 6 near midnight and 7 near noon 
8 and 9 were also similar, the two flights occurring at approximately the same local time, two days apart. 
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Figure 22 Calculated temperature profiles for an auroral arc with A3914 A brightness of 34 7 kR. 
[T 1 is the result neglecting thermal conduction, from Walker and Rees (1968a)] 
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Figure 23. Calculated temperature profiles tn a SACARC when the heating is supposed 
to result from (a) electric fields of various magnitudes normal to the magnetic field and 
(b) heat fluxes conducted from the magnetosphere. The temperature gradients are those 
at 1000 km. [From Walker and Rees (1968b)] 
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